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I  I.  INTRODUCTION 

^  This  document  describes  several  programs  used  by  MARC  in  studying  and  testing 
confidence  ellipses..  All  of  these  programs  were  developed  on  a  Hewlett 
Packard  Series  200  cmputer,  the  HP9836.  and  are  written  in  Hewlett  Packard's 

language.  This  a  very  powerful  version  of  BASIC,  and 
the  programming  constructs  found  in  these  programs  are 
non-standar(i«  Also,  the  graphics  commands  are  peculiar  to  this  BASIC. 


BASIC  3>0  prograimnii 
consequently  man)^<Jf 
U  uso. 


and  Ellipstuff.  The  listings  and 
The  programs  themselves  are 


pi^'^^ograms  are  called  Genelllpse,  Ellsim, 
comments  are  current  as  of  August  9,  1985. 
subject  to  change  without  notice. 

^Genelllpse  is  an  "ellipse  graphing  and  combination'"  program.  It  allows  the 
user  to  specify  up  to  twenty  confidence  ellipses,  combine  ellipses,  graph 
ellipses,  and  test  ellipses  for  combination.  It  has  been  used  primarily  to 
explore  the  geometric  properties  of  the  combination  method,  but  has  also  been 
used  to  create  figures  for  a  number  of  other  ^ports. 


rm  'nr 


Ellsim  is  a  "confidence  ellipse  simulation"  program.  It  has  been  used  to 
explore  the  robustness  and  properties  of  the  statistical  test  (vdiich  is  used 
to  decide  whether  or  not  to  combine  two  ellipses).  In  it,  two  normal  data 
distributions  are  specified,  corresponding  to  two  emitters.  Confidence 
ellipses  are  generated,  tested,  and  combined,  and  various  descriptive 
statistics  are  compiled.  In  addition,  there  is  a  routine  to  calculate  the 
power  of  the  statistical  test  in  certain  cases. 


9^ 

Ellipstuff  is  an  "ellipse  routine  library, 
routines  to  make  working  with  ellipses  easier.  As  such,  both  Genelllpse  and 
Ellsim  contain  all  of  these  routines,  but  they  are  also  listed  separately.  If 
more  programs  need  to  be  written,  these  are  the  routines  to  build  them  with. 


which  contains  a  large  number  of 


Each  of  these  programs  is  listed  below,  along  with  a  short  description  of 
their  use.  This  is  a  supplement  to  the  report  "Testing  and  Combination  of 
Confidence  Ellipses:  A  Geometric  Analysis,"  submitted  to  JPL  by  MARC  on 
August  5,  1985.  Refer  here  any  questions  concerning  confidence  ellipses,  the 
combination  method,  or  the  statistical  test  used. 


II 


GUIDE  TO  GENELLIPSE;  Ellipse  Graphing  Program 


Genellipse  stands  for  General  Ellipse  graphing,  and  can  also  combine  ellipses 
and  test  two  ellipses  for  combination.  The  program  makes  extensive  use  of  the 
Ellipstuff  Library  Routines. 

VRien  Genellipse  is  executed,  the  following  menu  is  displayed  on  the  screen: 


GENERAL  ELLIPSE  GRAPHER 


(A)  —  Enter  Ellipse  Mean  Point  (Center) 

(B)  —  Enter  Ellipse  Shape  (Covariance  Matrix) 

(C)  —  Enter  Ellipse  Shape  (Axes  and  Orientation) 

(D)  —  Combine  Two  Ellipse  (JPL  Method) 

(E)  —  Graph  an  Ellipse 

(F)  —  Clear  Graphics  Screen 

(G)  —  Choose  Plotter 

(H)  —  Set  Graphics  Screen  Bounds 

(I)  —  Calculate  Acceptance  Test  Statistic  (Chl**square) 

(J)  —  Draw  Axes 
(X)  —  Exit  Program 

Enter  your  choice: 


These  options  will  be  dealt  with  one  by  one. 

(A)  —  Enter  Ellipse  Mean  Point  (Center) 

This  option  prompts  the  user  to  enter  the  center  point  for  one  or  more 
ellipses.  The  program  first  asks  for  the  ellipse  (1-20),  and  then  for 
the  mean  point.  This  will  repeat  until  all  means  have  been  entered. 
When  the  user  is  through  entering  means,  pressing  the  return  key  in 
response  to  the  question  "Getting  mean  point  for  Ellipse  #"  will  return 
the  user  to  the  menu  shown  above. 


(B)  —  Enter  Ellipse  Shape  (Covariance  Matrix) 

There  are  two  ways  of  specifying  the  shape  of  a  confidence  ellipse.  The 
first  is  through  a  covariance  matrix.  This  option  will  ask  for  the 
ellipse  number,  and  then  prompt  for  the  elements  of  the  matrix.  It 
repeats  until  all  matrices  have  been  entered.  Just  as  option  A  does. 


(C)  —  Enter  Ellipse  Shape  (Axes  and  Orientation) 

The  other  way  of  specifying  the  shape  of  a  confidence  ellipse  Is 
geometrically.  This  option  prompts  for  the  ellipse  number.  Just  as  those 
above,  and  then  asks  for  the  axes  lengths  and  orientation  in  degrees  of 
the  ellipse.  It  repeats  until  all  ellipses  have  been  entered.  Just  as 
option  A  does. 

(D)  —  Combine  Two  Ellipses  (JPL  Method) 

This  option  combines  confidence  ellipses  using  the  method  described  in 
"Testing  and  Combination  of  Confidence  Ellipses:  A  Geometric  Analysis." 
The  option  prompts  for  the  numbers  of  the  two  ellipses  to  be  combined, 
and  then  for  the  number  of  the  resultant  ellipse.  For  example,  suppose 
the  user  had  entered  ellipses  1,  2,  3*  and  7  out  of  20,  and  wanted  to 
combine  2  and  3<  The  resultant  ellipse  could  be  given  number  <l,  5,  6,  or 
8  through  20.  Further,  if  either  of  ellipses  1  or  7  were  no  longer 
needed,  the  resultant  ellipse  could  be  given  either  of  numbers  1  or  7  as 
well.  This  would,  of  course,  erase  the  ellipses  originally  stored  in 
these  slots. 

(E)  —  Graph  an  Ellipse 

This  option  prompts  for  the  number  of  the  ellipse  the  user  wishes  to 
graph,  and  continues  asking  for  ellipses  until  all  desired  have  been 
graphed.  The  process  is  similar  to  that  in  option  A.  The  ellipse  will 
be  graphed  on  the  current  graphics  device  (See  option  G).  Other  options 
related  to  graphing  are  G,  H,  and  J. 

(F)  —  Clear  Graphics  Screen 

If  the  CRT  screen  is  the  current  graphics  device,  choosing  this  option 
will  erase  all  ellipses  currently  drawn  on  it. 

(G)  —  Choose  Plotter 

This  options  allows  the  user  to  select  the  current  graphics  device.  This 
program  is  currently  written  to  graph  on  the  CRT  screen  and  on  an  HP7470A 
two-pen  plotter.  Consequently,  if  selected,  the  program  will  ask  "Do  you 
want  to  (D)raw,  (E)rase,  Plot  with  Pen  (0),  (1),  (2)?"  Choosing  "D"  will 
cause  it  to  do  all  graphing  on  the  screen.  Choosing  "E"  will  cause  it  to 
graph  on  the  screen  in  "black"  —  so  that  an  ellipse  may  be  erased 
without  clearing  the  whole  screen.  Choosing  (1)  or  (2)  will  cause  it  to 
graph  on  the  plotter,  using  the  specified  pen.  Choosing  (0)  causes  the 
plotter  to  put  away  the  pen  that  its  using,  and  then  sets  the  CRT  screen 
to  be  the  graphics  device. 

(H)  —  Set  Graphics  Screen  Bounds 

This  option  is  used  to  scale  the  graphics  screen.  Default  scaling  is  -50 
to  50  on  the  X-axis  and  -50  to  50  on  the  Y-axis.  When  this  option  is 
chosen,  it  will  first  present  the  "limits"  of  any  ellipse:  that  is,  how 
far  the  ellipse  extends  in  the  X  and  Y  directions,  so  that  a  reasonable 
screen  size  may  be  chosen.  The  program  prompts  for  the  number  of  each 


ellipse.  Just  as  in  option  A.  Press  the  return  key  after  all  desired 
ellipse  limits  have  been  seen.  The  program  will  then  ask  for  the  screen 
bounds:  Minimum  X  value,  Maximum  X  value,  Minimum  Y  value.  Maximum  Y 
value. 

-  Calculate  Acceptance  Test  Statistic  (Chi-square) 

This  option  prompts  the  user  for  two  ellipse  numbers,  and  then  performs 
the  statistical  test  described  in  "Testing  and  Combining  Confidence 
Ellipses:  A  Geometrical  Analysis."  The  value  of  the  test  statistic  is 
printed  on  the  screen.  If  it  is  less  than  or  equal  to  the  95}  Chi-square 
value,  5.991,  then  the  test  accepts;  otherwise  the  test  rejects. 
(Actually,  the  test  may  be  run  at  any  confidence  level  desired.  However, 
the  ellipses  used  in  this  program  are  assumed  to  be  95}  confidence 
ellipses,  and  if  the  test  is  performed  at  any  other  level  the  geometrical 
results  shown  in  "Testing  and  Combining..."  concerning  the  test  will  not 
necessarily  hold  true. 

-  Draw  Axes 

This  option  draws  a  set  of  axes  on  the  current  graphics  device,  and 
labels  them  according  to  the  current  screen  bounds  (see  Options  G  and  H). 


III.  GUIDE  TO  ELLSIM;  Ellipse  Simulation  Program 


Ellslm  stands  for  Ellipse  Simulator,  but  can  also  find  analytical  estimates 
of  statistical  power.  It  uses  the  Ellipstuff  Library  extensively. 

See  the  report  "Testing  and  Combining  Confidence  Ellipses:  A  Geometric 
Analysis"  for  a  description  of  simulation  and  how  it  has  been  applied  to  the 
ellipse  combination  problem.  This  is  the  program  used  to  generate  the 
results  in  section  V  of  that  report.  Mote  that  in  this  program,  sensor  error 
is  assumed  to  follow  a  bivariate  normal  distribution  about  the  true  location 
of  the  emitter.  When  location  estimates  are  derived  from  lines  of  bearing, 
however,  this  assumption  may  be  unrealistic. 

When  the  program  is  run,  it  will  present  the  user  with  the  following  menu  of 
choices: 


Ellipse  Combination  Program  Driver 

(A)  —  Specify  True  Covariance  Matrices 

(B)  —  Specify  True  Mean  Parameters 

(C)  —  Specify  Observations 

(D)  —  Call  Simulation  Generator 

(E)  —  Call  Power  Generator 
(X)  —  Exit  Program 


Enter  your  choice: 


Options  A,  B,  and  C  are  used  to  specify  the  two  data  distributions. 

Option  A  prompts  for  the  means  of  two  bivariate  normal  distributions. 

Setting  the  means  to  be  equal  is  equivalent  to  having  only  one  emitter. 
Setting  them  apart  is  equivalent  to  having  two  emitters. 

Option  B  prompts  for  the  covariance  matrices  for  the  two  bivariate  normal 
distributions. 

Option  C  prompts  for  the  sample  sizes  to  be  used  for  each  distribution.  Note 
that  any  confidence  ellipses  generated  will  have  covariance  matrices  equal  to 
those  specified  in  option  B,  divided  by  these  sample  sizes. 

Options  D  and  E  call  the  Simulation  Generator  and  the  Power  Generator 
respectively.  These  will  be  discussed  individually. 


(D)  —  The  Simulation  Generator 


Choosing  this  option  leads  to  another  menu:  the  simulation  menu.  It  appears 
as  follows: 


Enter  the  Letter  of  your  choice: 

(A)  Reset  the  Random  Number  Seed 

(B)  Simulate  using  true  Covariance  Matrices 

(C)  Simulate  using  estimated  Covariance  Matrices 

(D)  Display  Results  on  the  Printer 
(X)  Exit  Program 


Enter  Your  Choice: 


These  options  will  be  dealt  with  one  by  one 


(A)  —  Reset  the  Random  Number  Seed 

This  option  allows  the  user  to  start  the  random  number  process  with  a 
given  seed;  this  is  generally  not  necessary,  but  is  useful  for  purposes 
of  debugging. 

(B)  —  Simulate  using  true  Covariance  Matrices 

This  option  will  prompt  the  user  for  the  number  of  simulations  (generally 
100  or  more),  and  then  proceeds  in  this  manner:  for  each  simulation  it 
generates  the  number  of  observations  specified  by  the  sample  size  entered 
earlier.  It  estimates  the  emitter  locations  from  these  observations,  and 
calculates  confidence  ellipses.  It  tests  to  see  if  these  ellipses  may  be 
combined.  Finally,  it  checks  to  see  if  the  combined  ellipse  contains  the 
true  location  of  the  emitterCs).  When  all  simulations  are  done,  it 
compiles  these  results.  See  "Testing  and  Combining  Confidence  Ellipses" 
for  more  information. 

(C)  —  Simulate  using  estimated  Covariance  Matrices 

The  process  here  Is  the  same  as  that  outlined  for  option  B,  except  for 
one  addition:  the  covariance  matrices  used  in  the  confidence  ellipses 
and  statistical  test  are  estimated.  However,  the  formulas  for  the 
ellipses  and  test  assume  that  the  covariance  matrices  are  known.  Thus, 
this  option  is  used  to  explore  what  happens  if  estimated  matrices  are 

mistakenly  used.  Note  that  the  estimates  are  made  using  the  "S^” 
statistic.  This  is  the  usual  way  of  estimating  variance-covariance  from 
a  data  set,  but  is  different  from  the  methods  used  in  most  if  not  all  of 
the  position  fixing  algorithms  we  have  seen. 


ikiMilSLhM 


(D)  —  Display  Results  on  the  Printer 

When  either  option  B  or  option  C  has  been  completed,  the  results  of  the 
simulations  are  shown  In  the  screen.  If  a  hardcopy  Is  desired,  selecting 
this  option  will  cause  the  results  of  the  last  simulation  run  to  be 
output  on  the  printer. 

(X)  —  Exit  Program 

This  option  will  exit  the  Simulation  Generator,  and  return  to  the 
original  menu. 
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(E)  The  Power  Generator 


The  power  of  a  statistical  teat  is  essentially  the  probability  that  the  test 
will  reject  when  it  ought  to.  That  is,  it  is  the  probability  that  the 
statistical  test  will  say  that  there  are  two  emitters  when  in  fact  there  are 
two  emitters.  In  the  problem  at  hand,  however,  the  power  is  not  a  single 
quantity:  in  fact,  there  is  a  different  power  value  for  each  pair  of 
emitters.  If  the  emitters  are  close  together,  the  power  of  the  test  will  be 
low;  if  they  are  far  apart,  the  power  will  be  close  to  1 . 

This  option  works  in  the  following  way,  for  convenience  sake.  It  uses  the 
covariance  matrices  and  sample  sizes  specified  from  the  main  menu,  but  allows 
the  user  to  enter  the  distance  between  the  two  emitters.  It  calculates  the 
power,  and  then  asks  for  another  distance.  To  return  to  the  main  menu,  enter 
0  for  the  distance. 


IV.  GUIDE  TO  ELLIPSTUFF;  Ellipse  Routine  Library 


Ellipstuff  is  a  library  of  subprograms  and  functions  written  to  facilitate 
confidence  ellipse  research  for  JPL.  It  covers  such  things  as  defining 
ellipses,  combining  ellipses,  testing  ellipses  for  combination,  and  graphing 
ellipses. 

To  use  Ellipstuff,  simply  Include  it  in  within  a  progreun.  Genellipse  and 
Ellslm  are  examples  of  this. 


Confidence  ellipses  and  how  they  are  stored 

A  confidence  ellipse  is  defined  by  two  things:  a  point  estimate,  or  mean, 
and  a  covariance  matrix.  Thus,  both  of  these  pieces  of  information  must  be 
stored  for  each  ellipse.  In  addition,  it  is  often  necessary  to  have  the 
inverse  of  the  covariance  matrix  on  hand  as  well.  Ellipstuff  stores  ellipses 
in  a  matrix  with  40  slots,  allowing  the  storage  of  20  covariance  matrices 
with  their  Inverses.  In  general  most  of  the  Ellipstuff  routines  deal  with 
the  inverses  themselves,  but  in  case  it  is  necessary  to  use  them  explicitly, 
the  ellipses  are  stored  in  slots  1  through  20  and  the  inverse  covariances  are 
stored  in  slots  21  through  40.  By  convention,  the  inverse  of  the  covariance 
matrix  for  the  ellipse  in,  say,  slot  3,  is  stored  in  slot  23 •  and  so  on. 
Examine  the  routines  for  more  programming  information. 


Ellipstuff  User  Routines 

The  slots  specified  in  the  following  routines  should  be  between  1  and  20 
Inclusive. 

Get_ell  jDean( Ellipse ) 

This  command  prompts  the  user  to  input  the  mean  (center  point)  for 
the  ellipse  in  slot  Ellipse. 

Example:  Getjelljneand ) 

Getjsovar lance ( Ellipse ) 

This  command  prompts  the  user  to  input  the  covariance  matrix  for  the 
ellipse  in  slot  Ellipse. 

Example:  Getjeovarlance(2) 

Get  jixes ( Ellipse ) 

This  command  prompts  the  user  to  input  the  shape  of  the  ellipse 
(lengths  of  the  semi -minor  and  semi-major  axes,  and  the  orientation) 
tdiich  is  then  converted  to  a  covariance  matrix. 

Example:  Get  axes(3) 


Test(Elllpse1 ,  Elllpse2,  Work,  Test_stat} 

This  command  runs  a  chi-square  test  on  the  ellipses  in  slots 
Ellipsel  and  EllipseE  in  order  to  see  if  they  may  be  combined.  The 
value  of  the  test-statistic  is  returned  in  Test__stat.  Work  is  the 
number  of  any  unused  slot,  to  be  used  for  scratch  work. 

Example:  Test(3t4t20,Some_variable) 

Combinejelllpse(Elllpae1 ,  EllipseE,  Combo) 

This  command  combines  the  ellipses  and  point  estimates  in  slots 
Ellipsel  and  Ellipse2,  and  stores  the  combined  ellipse  in  slot 
Combo. 

Example:  Comblne_ellipse(3,4,5) 

Draw_ellipse(Elllpse,  Xmin,  Xmax,  Ymin,  Ymax.Probjsonstant) 

This  draws  the  ellipse  in  slot  Ellipse  on  the  current  plotting 
device  (see  Choose_j>lotter).  Xmin,  Xmax,  Ymin,  and  Ymax  specify  the 
screen  dimensions.  The  default  values  are  (-50,50,-50,50). 
Prob_constant  is  zero  minus  the  chi-square  cutoff  associated  with 
the  confidence  level.  (The  cutoff  is  P--2LN(1 -Alpha) ,  where  Alpha 
is  the  confidence  level.  For  some  reason,  all  of  the  programming 
was  done  in  terms  of  -P-2LN(1 -Alpha).  Thus,  for  95K  ellipses, 
Prob_constant  -  -5.991). 

Example:  Draw__elllpse(2,-10,10,-20,20,-5.991 ) 

Disp_extremes(Elllpse) 

This  command  displays  the  extreme  x  and  y  points  calculated  by 
Getjoound  for  the  ellipse  in  slot  Ellipse.  This  is  designed  to  aid 
in  choosing  the  Xmin,  Xmax,  Ymin,  and  Ymax  values  required  by 
Drawjsllipse. 

Exan^le:  DiSp__extremes(1 ) 

D i splay_cov ( Ell i p se ) 

This  command  displays  the  covariance  matrix  of  the  ellipse  in  slot 
Ellipse  on  the  crt  screen  (not  graphically).  It  may  also  be  used  to 
display  the  inverses  in  slots  21  through  <t0. 

Example :  Display_cov (11) 

Choose_plotter 

When  first  initialized,  Ellipstuff  assumes  that  all  graphing  will  be 
done  on  the  CRT.  Choose^lotter  is  called  to  allow  the  program  user 
to  select  iriiich  device  to  use.  Note:  this  program  was  written  on  a 
system  with  two  graphics  devices — a  CRT  screen  and  an  HP7470A  two 
pen  plotter.  Choose_plotter  presents  the  user  with  5  options: 

D) raw  —  graph  on  the  CRT  in  the  normal  fashion  (DEFAULT). 

E) rase  —  graph  on  the  CRT  in  "black".  This  may  be  used  to  erase 
things. 

1 )  graph  on  the  plotter  using  Pen  1 . 

2)  —  graph  on  the  plotter  using  Pen  2. 

0)  —  put  away  the  current  plotter  pen  and  select  the  CRT. 

Example:  Choose_plotter 


Elllpatuff  Low  Level  Routines 

These  are  routines  which  are  used  in  building  user  routines.  Since  these 
oust  deal  with  Inverse  covariance  matrices  also,  the  slots  specified  may  run 
from  1  to  40.  The  user  will  generally  not  use  these,  except  when  adding  to 
the  Ellipstuff  module. 

Invert (Source,  Destination) 

This  command  Inverts  the  covariance  matrix  in  slot  Source  and  puts 
the  Inverse  in  slot  Destination. 

Example:  InvertOt  23) 

• 

Add_covariance(Elllpse1 ,  Ellipse2,  Summer) 

This  command  adds  the  covariance  matrices  stored  in  slots  Ellipsel 
and  Elllpse2,  and  stores  them  in  slot  Summer. 

Example:  Add_covar lanced ,2,3) 

Get_bounds(Ellipae ,  Probjconstant) 

~  This  command  calculates  and  saves  the  extreme  x  and  y  points  for  the 
ellipse  in  slot  Ellipse. 

Example:  Get  bounds(3«  -5.991) 
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GENELLIPSE: 


Ellipse  Graphing  Program 


Genelllpse  Is  an  "ellipse  graphing  and  combination"  program.  It  allows  the 
user  to  specify  up  to  twenty  confidence  ellipses,  combine  ellipses,  graph 
ellipses,  and  test  ellipses  for  combination.  It  has  been  used  primarily  to 
explore  the  geometric  properties  of  the  combination  method,  but  has  also  been 
used  to  create  figures  for  a  number  of  other  reports. 


GENERAL  ELLIPSE  GRAPHER 


100  !•' 

110  !• 

120  l» 

130  t*  This  prograo  uses  the  general  ellipse  generation  and  graph-*! 

100  I*  Ing  routines  found  In  the  file  ELLIPSTUFF.  Its  purpose  ■! 

ISO  I*  is  to  graph  general  ellipses  and  combinations  of  ellipses.  *! 

160  !•  •! 

170  I*  Hill  Duquette  Hay  13.  198S  *1 

190  COM  /Ellipses/  Matrx(2.2.il0).Xbar(20).ybar(20),Hajor_axis(20).Mlnor_axis( 

20),Theta(20).Angle(20).Ealn_x(20).Emax_x(20),Enln_y(20).Ew_y(20)  ~ 

191  REAL  Xaln.Xnax.ynln.YnSu.Probjsenstant 

200  Inlt  prog:!  ~ 

210  gInIT 

220  GCLEAR 

230  !  Ihese  variables  define  the  size  of  the  screen.  For  beat  results, 

2<I0  !  the  lengths  of  the  X  and  X  axes  should  be  about  the  same. 

250  Xain«-50 

260  XBax>50 

270  Inin— so 

280  Xnax>50 

281  Probjsonatant>2*L0G( .05) 

290  I  oTsplay  the  menu 

300  Henu_top:! 

310  Clear screen 

320  PRINT  "GENERAL  ELLIPSE  GRAPHER" 

330  PRINT 

3«0  PRINT  — - — - " 

350  PRINT  "  (A)  —  Enter  Ellipse  Mean  Point  (Center)" 

360  PRINT  "  (B)  —  Enter  Ellipse  Shape  (Covariance  Matrix)" 

370  PRINT  "  (C)  —  Enter  Ellipse  Shape  (Axes  and  Orientation)" 

380  PRINT  "  (0)  —  CoAine  Two  Ellipses  (JK,  Method)" 

390  PRINT  "  (E)  —  Graph  an  Ellipse” 

too  PRINT  "  (F)  —  Clear  Graphics  Screen" 

*01  PRINT  "  (G)  —  Choose  Plotter" 

*10  PRINT  "  (H)  —  Set  Graphics  Screen  Bounds" 

*20  PRINT  "  (I)  —  Calculate  Acceptance  Test  Statistic  (Chi-square)” 

*21  PRINT  "  (J)  —  Draw  Axes" 

*30  PRINT  "  (X)  —  Exit  Program" 

**0  PRINT 

*50  PRINT  "Enter  your  oholoe:" 

*60  OISP  "Choose  an  option:" t 

*70  INPUT  Options 

*80  SELECT  Options 

*90  CASE  "A", "a"  1  Get  ellipse  mean  points. 

500  Clearsoreen 

510  PRINT  "Getting  mean  point  for  Ellipse#"; 

520  REPEAT 

530  Elllpae«0 

5*0  DISP  "Enter  the  Ellipse#  (1-20)"; 

550  INPUT  Ellipse 

560  PRINT  Ellipse 

570  IF  Elllpae<-20  AMO  Elllpse>0  THEN  CALL  Get  ell  nean(Elllpae) 

580  UNTIL  Elllpse-0  ~  ~ 

590  CASE  ••B"."b”  I  Get  ellipse  covariance  aatricea 

600  Clearsoreen 

610  PRINT  "Getting  Shape  (Covariance  Matrix)  for  Ellipse#"; 

620  REPUT 

630  Ellipse-0 

6*0  DISP  "Enter  the  Ellipse#  (1-20)"; 

650  INPUT  Ullpse 

660  PRINT  Ullpse 

670  IF  Ullpae<-20  AMD  Ullpse>0  THEN  CALL  Get  cover lanee(Elllpae) 

680  UNTIL  Elllpse-0 

690  CASE  "C".”o"  I  Get  ellipse  axes  and  orientation 


700  ClMrscr««n 


710  PRINT  "C«ttlng  Shap«  (Axm  mui  Orientation)  for  EIllpMf": 

720  REPEAT 

730  ElllpM-0 

7«0  OISP  "Enter  the  Ellipse#  (1*20)*t 

750  INPUT  Ellipse 

760  PRINT  Ellipse 

770  IF  Elllpse<>20  AND  Elllpae>0  THEN  CALL  Get  axes(Elllpse) 

780  UNTIL  Elllpse-O 

790  CASE  1  Combine  ellipses 

600  Cleerscreen 

810  PRINT  "Combining  two  ellipses" 

820  PRINT 

830  Elllpse1-0 

640  Elllpse2-0 

850  Elllpse3-0 

860  PRINT  "Ellipse  1  #*: 

870  DISP  "Enter  the  first  Ellipse#  (1-20)"; 

880  INPUT  Elllpsel 

890  PRINT  Elllpsel 

900  PRINT  "ellipse  2  #"; 

910  OISP  "Enter  the  second  Ellipse#  (1-20)"; 

920  INPUT  Elllpse2 

930  PRINT  ElllpseE 

940  PRINT  "Combined  Ellipse  #*; 

950  OISP  "Enter  the  combined  Ellipse#  (1-20)"; 

960  INPUT  Ellipses 

970  PRINT  Ellipses 

980  IF  Elllpsel <-20  AND  Elllpsel >0  AND  Clllpse2<-20  AND  Elllpsel >0  AND  Ell 

lpse3<-20  AND  Clllpae3>0  AND  ElllpseSOElllpsel  AND  Elllpse3<>Elllpse2  THEN 
990  Combine  elllpse(Elllpse1 .Elllpse2, Ellipses) 

1000  END  IF  ~ 

1010  PRINT  "New  Heen  Point:  ",Xbar(Elllpse3).">  ",»ar(Elllpse3) 

1020  PRINT 

1030  PRINT  "New  Covariance  Matrix:" 

1040  Dlsplair_cov(Elllpse3) 

1050  PauseabTt 

1060  CASE  "E"."e"  I  Graph  ellipses 

1061  REPEAT 

1070  Clearscreen 

1080  PRINT  "Graphing  Ellipse#"; 

1100  Elllpae-0 

1110  OISP  "Enter  Ellipse#  (1-20)»; 

1120  INPUT  Ellipse 

1130  PRINT  Ellipse 

1140  IF  Elllpse<-20  AND  Elllpse>0  THEN 

1150  PRINT 

1160  PRINT  "Center  Point  («|Xbar(Elllpse) ;",  ";Ibar(Elllp8e) ;»)" 

1170  PRINT 

1180  PRINT  "Covariance  Matrix:* 

1190  Dlsplap  eov(Elllpse) 

1200  PRINT  "Major:  ";Major  axls(Elllpse) ;”  Minor:  *;Hlnor  axla(Elllp 

se)  ” 

1210  PRINT  "Orient:  ";Angle(Elllpse) 

1220  Invert(Elllpse,Elllpae*20) 

1230  Oee_bounds(Elllpse,Prob  constant) 

1250  Draw  ell Ipse (Ell ipse, XmTn.Xmax.Tmin.Xnax.Prob  oonstant) 

1260  ALPHA  OH 

1280  END  IF 

1281  UNTIL  Elllpse-O 

1290  CASE  ”F","r"  I  Clear  the  graphics  screen 

1300  QCLEAR 

1301  CASE  «C",*s"  I  Choose  the  plotter  device 

1302  Choose_plotter 

1310  CASE  *H","h"  I  Set  the  screen  boundaries 
1320  Clearscreen 
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1330  PRINT  "Setting  Screen  Boundaries....* 

13lt0  PRINT 

13SO  REPEAT 

1360  Elllpae-0 

1370  OISP  *Dlsplay  Extreme  Points  for  Ellipse*  (1-20)”; 

13SO  INPUT  Ellipse 

1390  IP  Elllpae>0  AND  Elllpse<-20  THEN 

IHOO  Invert(Elllpse,Elllpse*20) 

1 910  Cet_bounda(Elllpse.Prob_conatant) 

1920  Dls?  extremes(Elllpae)  ~ 

1930  END  IF  ~ 

1990  UNTIL  Elllpse-0 

1950  PRINT 

1960  Get  X:! 

1970  “  PRINT  "Minimum  X  value: 

1980  OISP  "Enter  the  minimum  X  value"; 

1990  INPUT  Xmln 

1500  PRINT  Xmln 

1510  PRINT  "Maximum  X  value: 

1520  OISP  "Enter  the  maximum  X  value”; 

1530  INPUT  Xmax 

1590  PRINT  Xmax 

1550  IF  Xmin>-Xmax  THEN  Get_x 

1560  Get_y:l  “ 

1570  PRINT  "Minimum  T  value:  "; 

1560  OISP  "Enter  the  minimum  T  value"; 

1590  INPUT  rmln 

1600  PRINT  Xmln 

1610  PRINT  "Maximum  X  value:  "; 

1620  OISP  "Enter  the  maximum  X  value"; 

1630  INPUT  Xmax 

1690  PRINT  Xmax 

1650  IF  Xmin>-Xmax  THEN  Get 

I660  CASE  "l","i"  1  Calculate  the  chi-square  acceptance  test  statistic 

1670  Clearsoreen 

1680  PRINT  "Calculating  Acceptance  Test" 

1690  PRINT 

1700  Elllpsel-O 

1710  ElllpseE-O 

1720  Swork-0 

1730  PRINT  "Ellipse  1 

1790  OISP  "Enter  the  first  Ellipse#  (1-20)"; 

1750  INPUT  Elllpael 

1760  PRINT  Elllpael 

1770  PRINT  "ellipse  2  *"; 

1780  OISP  "Enter  the  second  Ellipse*  (1-20)"; 

1790  INPUT  Elllpae2 

1800  PRINT  Elllpseg 

1810  PRINT  "Scratch  Work  *"; 

1820  OISP  "Enter  the  scratch  work  #  (1-20)"; 

1830  INPUT  Swork 

1890  PRINT  SMork 

1850  IP  Elllpaai<-20  AND  Elllpse1>0  AMD  ElllPae2<-20  AND  Elllpael >0  AND  Swo 
rk<-20  AND  Swork>0  AND  SworkOEllipsel  AND  SworkOElllpseE  THEN 
i860  Test(Elllpae1 ,Elllpse2,SHork>Test  stat) 

1870  END  IF  “ 

1880  PRINT 

1890  PRINT  "The  Test  statistic  is  ”:Tcst  stat 
1900  PRINT  " 

1910  PRINT 

1920  Pauaeablt 

1921  CASE  I  Draw  in  the  axes 

1922  GRAPHICS  ON 

1923  CSIZE  2 

1925  MOVE  0,0 

1926  DRAW  0.100 
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1928  DRAW  0,0 

1929  UBEL 

1930  LORC  7 

1931  HOVE  100,0 

1932  UBEL  Xmax 

1933  LORC  3 

1930  HOVE  0,100 

1935  UBEL  Xmax 

1936  LORC  V 

1937  HOVE  50,0 

1938  CSIZE  3 

19X0  UBEL  "X-Axia  (in  kiloaetera)* 

1941  HOVE  0,50 

1942  DEC 

1943  LDIR  270 

1944  UBEL  "r-Axia  (in  klloMtera)* 

1945  LORC  1 

1946  LOIR  0 

1947  RAD 

1949  CASE  "X","x"  I  Wa  can  atop  nou 

1950  Claaracreon 

1951  PRINT  "Thaf  a  all,  folkal" 

1960  STOP 

1970  CASE  ELSE 
1980  PRINT  CHR$(7) 

1990  END  SELECT 
2000  COTO  Hanu  top 
2010  END  “ 

5000  1 

5010  i  SUBROUTINES:  Takan  from  ELLtPSTUFF 

5020  I 

5030  SUB  In»art(Srca,Daat) 

5040  t  Thla  routlna  Invarta  any  eovarlanea  natrlx  In  Hatrx  and  placaa 
5050  I  tha  invartad  ouitrlx  In  Daat. 

5060  COH  /Clllpaaa/  Hatrx(2,2,40),Xbar{20),Ybar(20).HaJor  axla(20),Hlnor  ax 

la  ( 20 ) ,  Tha  ta  ( 20 ) ,  Angla(  20 ) ,  Enln_x(  20 ) ,  Eisax jt{  20 ) ,  Enln_y  ( 20 ) ,  Eamx_y  ( 20 ) 

5070  0at>Hatrx(1 ,1 ,Srca)*Natrx(2,2,Srca)>Hatrx(1 ,2,Srca)*HaTrx( 1 ,2,Srca) 

5080  Katrxd ,1,Deat)>Hatrx(2,2,Srca}/Dat 

5090  Hatrx(2.2,Da8t)>Hatrx(1 ,1  ,Srca)/Dat 

5100  Hatrxd ,2,0aat)«*Hatrx(1 ,2,Srca)/Dat 

5110  Hatrx(2,1 ,  Daat  )»'-Hatrx(  2,1  ,Srca)/Dat 

5120  SUBEND  I  End  of  SUB  Invar t 

5130  I 

5140  I  CET  BOUNDS 
5150  I 

5160  SUB  Cat  baunda(Elllpaa.Prob  oonstant) 

5170  I  iRla  aubroutlna  ealouXatas  tha  X  and  T  llalta  for  tha  glvan  alXlpaa 

5180  COH  /tlllpaaa/  Natrx(2,2,40),Xbar(20).Xbar(20),Kajor  axla(20).Hlnor  ax 

la(20),Thata(20),An«la(20),EBin_x(20),Eaax_x(20),Baln y(20),taax  t(20)  ~ 

5190  REAL  Taap1,Toap2,ToMpT.Taap4,Taiip5 

5210  Taapl*4latrx(2,2,Elllpao*20)*Prob  oonatant 

5220  Taap2«Hatrx(  1 ,2,Blllpaa*20)*Natrxd  ,2.Xlllpaa4’20)-Hatrx(  1 , 1  ,Elllpaa*20 

)*Hatrx(2,2,Elllpsa*20) 

5230  Taiip3-(Toap1/Ta«i2)‘(«5) 

5240  Toap4>Hatrx(1 ,1 ,Elllpaa*20)*Prob  oonatant 

5250  Taa95-(Tov4/Ta«2)*(.5) 

5260  lklnjr( Ell ipaa ) •-Toflp5*Tbar ( Ell ipao ) 

5270  Eaax_p(Ellipao)>Taap5*Vbar(Elllp8a) 

5280  lain ji( Ell Ipaa ) — Taap3*Xbar ( Ell ipaa  > 

5290  BBBx'‘x(Elllpaa)>Toap3*Xbar(Xlllpsa) 

5300  SUBEND  ~l  End  Of  GET  BOUNDS 
5310  I 

5320  I  CHOOSE  PLOTTER 


53*>0  SUB  Cl>ooa0_pl otter 

5350  I  Subroutine  to  ehooae  the  tfealred  plotter  device 

5360  REPEAT 

5370  Gojjn—1 

5380  DISP  "Do  you  want  to  (Dlraw,  (Elraae,  Plot  with  Pen  (0),  (1),  (2)7" 

5390  LINPUT  MeaaageS 

5<t00  SELECT  MeaaageS 

5*110  CASE  «D","d" 

5*120  PLOTTER  IS  3. "INTERNAL" 

5*130  GRAPHICS  OH 

5*1*10  PEN  1 

5*150  Co  on»1 

5*151  CASE  “O" 

5952  PLOTTER  IS  705,"HPGL" 

5953  GRAPHICS  ON 

5959  PEN  0 

5955  Go  on>1 

5960  CASE  ""l" 

5970  PLOTTER  IS  705."HPGL" 

5980  GRAPHICS  ON 

5990  PEN  1 

5500  Go  on>1 

5501  CASE  ^2" 

5502  PLOTTER  IS  705."HPGL" 

5503  GRAPHICS  ON 

5509  PEN  2 

5505  Go  on>1 

5510  CASE  "E","e" 

5520  PLOTTER  IS  3, "INTERNAL" 

5530  GRAPHICS  ON 

5590  PEN  -1 

5550  Co  on»l 

5560  CASE  ELSE 

5570  PRINT  CHR*(7) 

5580  END  SELECT 

5590  UNTIL  Co  on-t 

5600  SUBEMD  ~ 

5610  I 

5620  I  GET  ELL  MEAN 
5630  J 

5690  SUB  Getjsll_Bean(Clllpae) 

5650  I  ThXa  aiTbroutlnea  prompta  the  uaer  for  the  laean  of  an  elllpae. 

5660  COH  /Elllpaea/  Matrx(2,2,90),Xbar(20).Tbar(20),HaJor_axla(20),Mlnor_ax 

l8(20),Theta(20),Angle(20),Eain_x(20),EBax_x(20),Ealn_jr(20>,Eaax_y(20)  ~ 

5670  Clearaereen  —  —  —  — 

5680  PRINT  TABXT( 1.10} ("WHAT  IS  X>BAR  for  ELLIPSE  ";Elllpae;"7"{ 

5690  INPUT  Xbar( Elllpae) 

5700  PRINT  Xbar(Elllp8e) 

5710  PRINT  TABXXC  1,11) ("WHAT  IS  T-BAR  for  ELLIPSE  "(Elllpae ("7"; 

5720  INPUT  Xbar(Elllpae) 

5730  PRINT  Tbar(Elllpae) 

5790  Pauaeablt 

5750  SUBOrO 
5760  I 

5770  I  GET  COVARIANCE 
5780  1 

5790  SUB  Cetjeovarlance(Elllp8a) 

5800  I  liTla  routine  geca  the  covariance  natrlx  for  an  elllpae 

5810  CCH  /Ellipses/  Hatrx(2,2,90).Xbar(20).Xbar(20),HaJor_axla(20),Hlnor_ax 

la(20)  ,Theta(20)  ,Ansle(20)  ,BBln_x(20),EBax_x(20)  ,Eain_jr<20),Eaax_y(20) 

5820  Clearaereen  ~  ~ 

5830  FOR  K-1  TO  2 

5890  PRINT  TABXT(1,3*K*6)i"ENTER  ELEMENT  ("(Ki","{K(")  IN  THE  COVARIANCE 

MATRIX  FOR  ELLIPSEf"{Elllpse(":  "( 

5850  INPUT  Matrx(K.K,E11108e) 
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5860  PRINT  Matrx(K,K, Ellipse) 

5870  NEXT  K 

5880  PRINT  T*BXY(1,3*9);"ENTER  ELEMENT  (1,2)  IN  THE  COVARIANCE  MATRIX  FOR  E 
LEMENT#";!;"; 

5890  INPUT  Matrx( 1,2, Ellipse) 

5900  Hatrx(2,l ,Elllpse)-Matrx(1, 2, Ellipse) 

5910  PRINT  Hatrx( 1,2, Ellipse) 

5920  Pauseablt 

5930  SUBQID 
59*10  I 

5950  I  PAUSEABIT 
5960  i 

5970  SUB  Pauseablt 

5980  I  Pause  and  wait  for  a  carrlaRe  return 

5990  OISP  "Type  ENTER  to  continue..."; 

6000  INPUT  Garbaget 

6010  SUBEND 
6020  SUB  Clearscreen 
6030  I  Clear  the  screen 

60*10  PRINT  CHR$(12) 

6050  SUBEMD 
6060  I 

6070  I  GET  AXES 
6080  !  ~ 

6090  SUB  Get_axes(Elllpse) 

6100  I  This  subroutine  gets  an  ellipse  In  terms  of  the  axes  and  the 

6110  I  angle  of  orientation.  These  are  converted  into  a  covariance 

6120  i  matrix. 

6130  COM  /Ellipses/  Matrx(2,2,*l0),Xbar(20),Ybar(20),Major_axis(20),Mlnor_ax 

la(20),Theta(20),Angle(20),Emin_x(20),Emax  x(20),Emlnjr(20),Emax  y(20)  ~ 

61*10  Prob_constant«2*L0C(.05)  —  —  — 

6150  CleaFsereen 

6160  PRINT  "WHAT  IS  THE  RADIUS  OF  THE  MAJOR  AXIS  OF  ELLIPSE  ";Elllpse;"  s  " 

6170  INPUT  Major  axla(Elllpse) 

6180  PRINT  Major~axia( Ellipse) 

6190  PRINT  "WHAT~IS  THE  RADIUS  OF  THE  MINOR  AXIS  OF  ELLIPSE  "sElllpse;"  :  " 

6200  INPUT  Minor  axlsCEllipse) 

6210  PRINT  Hlnor'~axls(Elllpse) 

6220  PRINT  "HHAT~ IS  THE  ANGLE  (IN  DEGREES)  OF  THE  MAJOR  AXIS  OF  ELLIPSE  ";E 

lllpses"  s  "; 

6230  INPUT  Angle(Elllpse) 

62*)0  PRINT  Angle( Ellipse) 

6250  I  Convert  to  Covariance  Matrix 

6260  Theta(Ellipse)-PI"Angle(Elllpae)/180 

6270  Matrxd  ,1  ,Elllpse)-(Major  axia(Elllpse)*C0S(TheU(Elllpse)))*2*(Hlnor 

axls( Ellipse) *SIN(Theta(Elllpse)))*?  ~ 

6280  Natrx(2, 2, Ellipse) •(Major  axls(Elllp8e)*SIN(Theta(Elllpse)))*2»(Mlnor_ 

axis(Elllpse)*COS(Theta(Elllpse)))*?  ~ 

6290  Hatrx(2,1,Elllpse)«((MaJor  axls(Elllpse))*2*(Mlnor  axls(Elllpse) )*2)*C 

OS(TheU(Ellipse))*SIN(Theta(ElllpaeT)  ~ 

6300  Matrxd  ,2, Elllpse)>Hatrx(2,1  .Ellipse) 

6310  FOR  J-1  TO  2 

6320  FOR  K-1  TO  2 

6330  NatrxC J,K,Elllpse)>Matrx{J,K,£lllpse)/(-Prob  constant) 

63*10  NEXT  K 

6350  NEXT  J 

6360  PRINT 

6370  Pauseablt 

6380  SUBEMD 
6390  I 

6400  t  DISPUX  COV 
6410  I  ^ 

6420  SUB  Dlsplayj!Ov(Elllpse) 

6430  t  Displays  the  covariance  Mtrlx  for  an  ellipse 


6*tU0  COM  /Ellipses/  Matrx(2,2.<l0),Xbar(20),Xbar(20),Major_axls(20),Minor_ax 
la(20) .'nieta(20) .Angle(20)  .Emin  x(20),Emax  x(20),Emln  y(20),Emax  y(20) 

61*50  FOR  1-1  TO  2  ~  “  ~  “ 

6<i60  FOR  J-1  TO  2 

6^70  PRI.VT  Matrxd.J, Ellipse)!" 

6«80  NEXT  J 

6490  PRINT 

6500  NEXT  I 

6510  SUBEMD 
6520  I 

6530  I  ADD  COVARIANCE 
6540  ! 

6550  SUB  Add_covarlance(Flrst. Second. Sunmer) 

6560  I  nTls  subroutine  can  be  used  to  add  Suimer«Flrat»Second 

6570  CON  /Ellipses/  Hatrx(2.2.40) .Xbar(20) .Tbar(20) .Major_axls(20) .Mlnor_ax 

ls(20).Theta(20).Angle(20).Efflin  x(20),Emax_x(20).Enln  y(2D).Emax_y(20) 

6580  INTEGER  I.J  ~  ~ 

6590  FOR  1-1  TO  2 

6600  FOR  J-1  TO  2 

661 0  Hatrx( I .J .Summer )-Matrx( I .J .First ) *Matrx( 1 . J .Second ) 

6620  NEXT  J 

6630  NEXT  I 

6640  SUBEHD 
6650  I 

6660  t  COMBINE  ELLIPSE 
6670  I  ~ 

6680  SUB  Cofflblne_elllpse(Flrst. Second. Combo) 

6690  I  This  routine  finds  the  "JPL"  combination  of  the  First  and  Second 

6700  I  ellipses. 

6710  COM  /Ellipses/  Matrx<2.2.40) .Xbar(20) .Ybar(20) .Major_axls(20) .Mlnor_ax 

ls(20).Theu(20).Angle(20).£;fflln  x(20).Emax  x(20).Emln  y(20).Emax  y(20)  ** 

6720  REAL  Sxl .3x2 

6730  I  Calculate  the  new  covariance  matrix. 

6740  Invert(Flr8t.Flrst+20) 

6750  Invert(Second.Second*20) 

6760  Add_covarlanee(Flrst'»20.Second'»20,Combo*20) 

6770  Invert(Combo-20, Combo) 

6780  I  Calculate  the  new  mean  point. 

6790  Sx1-Matrx(1 .1 .Flrst*20)*Xbar(Flrst)*Matrx(1 .2,Flrst*20)»rbar(Flrst) 

6800  Sx1-Sx1*MatrxO  .1  .Second-20)*Xbar(Second)'»Matrx(1 .2,Second*20)*Ibar(Se 

cond) 

6810  Sx2-Matrx(2.1 .Flrst»20)»Xbar(Flrst)*Matrx(2.2,Flrst»20)»rbar(FIrst) 

6820  Sx2-Sx2*Natrx(  2 . 1 , Second-20 )*Xbar( Second ) -Matr x( 2 . 2 .Second-20 ) • Ybar  ( Se 

cond) 

6830  Xbar(Combo)-Hatrx(1 .1 .Combo)>Sxl-Matrx(1 .2.Combo)*Sx2 

6840  Ybar ( Combo ) -Natrx( 2 . 1 .Combo) "Sxl -Hatrx( 2.2. Combo) *Sx2 

6850  SUBOID 
6860  I 

6870  I  OISPJEXTREMES 
6880  I 

6890  SUB  01ap_oxtremea(Elllpse) 

6900  I  This  routine  displays  the  extreme  x  and  y  values  for  the  specified 
6910  I  ellipse 

6920  CON  /Ellipses/  Hatrx(2,2,40).Xt''?r(20),ybar(20).MaJor_axls(20).Mlnor_ax 

ls(20),Theta(20).Angle(20).Emln  x(20).Eaax_x(20).Emln  y(20).Emax  y(20) 

6930  PRINT  "E#  "iElllpse;"! 

6940  PRINT  "Xmln  ”;Emln  x( Ell Ipse); 

6950  PRINT  "!  Xmax  " iESix_x( Ellipse ) ; 

6960  PRINT  Tmln  ";Efflln  y( Ell ipse); 

6970  PRINT  tmax  "!Emax3<Elllpse) 

6980  SUBEMO 

6990  SUB  DrawjelllpseCElllpse. Xmln. Xmax. Ynln.Imax.Prob_constant) 

7000  I  ThlT  routine  draws  one  (1)  ellipse  on  the  current  plotter  device. 

7010  COM  /Ellipses/  Matrx(2.2.40) .Xbar(20) . Ybar(20) .Major_axls(20) .Mlnor_ax 

la(20).Theta(20).Angle(20).Emln  x(20).Emax_x(20).Emln  y(20) .Emax_y(20)  ~ 

7020  INTEGER  Sign 
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7030  REAL  Xpoint, Upolnt.Ypoint. Tempi ,Teop2, Temps, TeopU.TempS 

7060  GRAPHICS  ON 

7070  I  Draw  top  half  of  the  ellipse 

7080  Slgn-1 

7090  COSUB  Orawjialf 

7100  t  Draw  bottom  half  of  the  ellipse 

7110  Sign— 1 

7120  GOSUB  Draw_half 

7130  !  Finish  up 

71  *10  MOVE  ( Xbar  ( Ell  Ipse )  -Xmin )  •  1 00/  ( Xmax-Xmln ) .  ( Ybar  (Ellipse )  -Ynln ) » 1 00/  ( Ym 

ax-Ymln) 

7150  DRAW  ( Xbar ( Ellipse ) -Xmln ) *1 00/ ( Xmax-Xmln ) . ( Ybar ( Ell Ipse ) -Ymln ) • 1 00/ ( Ym 
ax-Ymln) 

7160  I  Okey,  dokey,  we're  done  here, 

7170  SUBEXIT 

7180  Draw_half:  I  Draw  half  of  the  ellipse 

7190  ~!  If  Slgn«l ,  then  draw  top  half;  If  Sign— 1 ,  then  bottom  half. 

7200  FOR  Xpolnt-Emln_x( Ellipse)  TO  Emax_x( Ell Ipse)  STEP  .1 

7210  (K)SUB  Draw_1_polnt 

7220  NEXT  Xpolnt 

7230  Xpolnt-Efflax_x( Ellipse) 

72*10  GOSUB  Draw_1_polnt 

7250  RETURN 

7260  1 

7270  1  Compute  each  point  and  draw  the  new  line.  (It's  here  since  we  call 

7280  1  It  twice 

7290  ! 

7300  Draw_1_point:  ! 

7310  Upolnt«Xpolnt-Xbar(Ellipse) 

7320  Tempi aMatrxd  ,2,Elllpae''20)*Upolnt 

7330  Temp2ti7emp1»Temp1-Matrx(2,2,Ellipse*20)*(Matrx(  1 ,1  ,Ellipse*20)»Upoint» 

Upolnt*Prob_constant) 

73*<0  IF  Temp2<10*(-10)  THEN  Temp2-0 

7350  Ypoint»(-Temp1+Slgn*SQR(Temp2))/Matrx(2,2,Elllpse*20) ♦Ybar  (Ellipse) 

7360  IF  Xpolnt<-Emin  x(Elllpse)».001  THEN 

7370  MOVE  (Xpoint~Xmin)»100/(Xmax-Xmln) , (Ypolnt-Ymin)»100/( Ymax-Ymin) 

7380  ELSE 

7390  DRAW  ( Xpoint-Xmin ) *1 00/ ( Xmax-Xmln ) , ( Ypoint-Ymln ) » 1 00/ ( Ymax-Ymin ) 

7*100  END  IF 

7*110  RETURN 

7*120  SUBEND 

7*130  I 

7*1*10  1  TEST 

7*150  ! 

7*160  SOB  Test(Firat, Second, Work, Test_val) 

7*170  I  This  routines  calculates  the  acceptance  test  criteria  for  First  and 

7*i80  I  Second.  Work  la  used  as  a  "scratchpad”. 

7*190  COM  /Ellipses/  Matrx(2,2,*IO) ,Xbar(20) , Ybar(20) ,MaJor_axis(20) ,Mlnor_ax 

ls(20),Theta(20),Angle(20),Efflln_x(20),Emax  x(20) ,Efflln_y(20) ,Emax_y(20) 

7500  Add_covariance( First, Second, WorkT 

7510  Invert(Nork,Work*20) 

7520  01ffx>Xbar(Flrst)-Xbar(Second) 

7530  Dlffy-Ybar(Flr8t)-Ybar(Seeond) 

75*10  Templ-Matrxd  ,1  ,Work^20)«Diffx*Matrx(1,2,Wopk^20)»Dlffy 

7550  Temp2«Matrx(2,1 ,WorK*20)«Dlffx*Matrx(2,2,Work*20)»Dlffy 

7560  Test_val-Dlffx*Templ ♦Dlffy"Temp2 

7570  SUBEND 
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ELLSIM: 


Ellipse  Simulation  Program 


Ellsim  is  a  "confidence  ellipse  simulation"  program.  It  has  been  used  to 
explore  the  robustness  and  properties  of  the  statistical  test  (which  is  used 
to  decide  whether  or  not  to  combine  two  ellipses).  In  it,  two  normal  data 
distributions  are  specified,  corresponding  to  two  emitters.  Confidence 
ellipses  are  generated,  tested,  and  combined,  and  various  descriptive 
statistics  are  compiled.  In  addition,  there  is  a  routine  to  calculate  the 
power  of  the  statistical  teat  in  certain  cases. 


1000 

1010 

loao 

1030 

10<I0 

1050 

1060 

1070 

1080 

1090 

1100 

1110 


ElllpM  Slaulatlon  Driver 

This  la  a  ample  almulatlon  program  for  the  elllpae 
combination  teat.  Only  the  Normal  dlatrlbutlon  la 
aupported. 

Original:  1/17/85  Updated:  7/29/85 


7/29/85 


! 


1120  Top  ofjirogram:  ! 

1130  ~l  Uae  Faat  Math  card 
1190  CONTROL  32.2 si 

1150  I  Specify  common  varlablea  between  the  different  aubprograma 

1160  COM  /Elllpaea/  Matrx(2,2,90),Xbar(20),ybar(20),MaJor_axla(20),Mlnor_axl 

a  ( 20 ) .  Theta  ( 20 ) ,  Angle(  20 ) ,  Emln_x(  20 ) ,  Emaxjt(  20 ) .  Emlnjr  ( 20 ) ,  Emaxjf  ( 20 ) 

1170  COM  /Driver/  0ba{2)  “ 

1180  Driver  menu:  i 


1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 


Clearacreen 

PRINT  "Elllpae  Combination  Program  Driver" 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


"(A)  —  Specify  True  Covariance  Matrlcea’ 
"(B)  “  Specify  True  ttean  Parametera" 

"(C)  —  Specify  Obaervatlona" 

"(D)  —  Call  Simulation  Generator" 

"(E)  —  Call  Power  Generator" 

"(X)  —  Exit  Program" 


1290  Get  choice: 
1300 
1310 


I 


INPUT  "Enter  your  choicej",0ptlon4 
SEUCT  Option* 

CASE  "A", "a" 

Clearacreen 

CALL  GetjJovarlanced  ) 

CALL  Cet'”covarianoe(2) 

CASE  "B"."F" 

Clearacreen 

CALL  Getjeil  mean(l) 

CALL  Get  elljBean(2) 

CASE  "C","o" 

Clearacreen 
GOSUB  Getjsba 
CASE  "0","d“ 

Clearacreen 
CALL  Sim  elllpae 
CASE  "E","e" 

Clearacreen 
CALL  Power  elllpae 
CASE  "X","x"~ 

GOTO  Endjirogram 
CASE  ELSE 

PRINT  CHH*{7) 

GOTO  Get  choice 
END  SELECT  ~ 

GOTO  Driver  menu 


1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 

1400 
1410 
1420 

1430 

1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 

1540 

1550 
1560 
1570 
1580 
1590 
1600 
1610 

1620  Get  oba:  I 
1630  “FOR  Imat-1  TO  2 


SUBROUTINES 


Get  the  number  of  TRIALS  (Main  Menu  Choice  "B") 


1640  PRINT  "How  oiany  observations  for  ellipse  #";Ioiat;": 

1650  INPUT  Obsdmat) 

1660  PRINT  Obsdmat) 

1670  NEXT  Imat 

1680  RETURN 
1690  End^program:! 

1700  END 

1720  I  SUBPROGRAMS  AND  FUNCTIONS 

1730  ! 

1740  SUB  Inver t(Srce,Deat) 

1750  1  This  routine  inverts  any  covariance  matrix  in  Hatrx  and  places 

1760  I  the  Inverted  matrix  in  Dest. 

1770  COM  /Ellipses/  Hatrx(2.2.40).Xbar(20).Ybar(20),HaJor_axis(20),Minor_ax 

la<20).Theta(20}.Angle(20).Ealn_x(20).Emax_x(20),Eain_y(20),Emax _/(20)  ~ 

1780  Det-Matrxd ,1 ,Sree)»Matrx(2,2,Sree)-Matrx(l72,Srce)»Matrx(1 .2,Srce) 

1790  Hatrxd  ,1 .0est)-Matrx(2,2.Srce)/Det 

1800  Hatrxd, 2. Best) •Hatrxd  .1  ,Sree)/Det 

1810  Hatrx( 1 .2,0est}>-Hatrx( 1 .2.Srce)/Det 

1 820  Hatrx( 2 . 1 . Dest ) — Hatrx( 2 . 1 ,Srce ) /Det 

1830  SUBEND  !  End  of  SUB  Invert 

1840  I 

1850  I  GET  BOUNDS 
1860  !  ~ 

1870  SUB  Get_bounds( Ellipse, Prob_conatant) 

1880  I  ItTis  subroutine  calcuTates  the  X  and  T  limits  for  the  given  ellipse 


1890  COM  /Ellipses/  Hatrx(2,2,4o).Xbar(20),Xbar(20),MaJor_axi8(20),Minor_ax 

ia(20),Theta(20),Angle(20),Emin  x(20),Emax_x(20),Eminjr(20) ,Emax_y(20)  ~ 

1900  REAL  Tempi  ,Temp2,TempT,TeHp4,Teii^5 

1910  Temp 1 •Matr x( 2,2, Ell ipse  *20 ) sprobjeonstant 

1920  Temp2«Matrxd  ,2,Ellipae*20)»Matr7d  ,2,Ellip8e*20)-Matrxd  ,1,Ellipse*20 

)  "Ma  trx(  2 , 2 ,  £11 1  paedO ) 

1930  Temp3-(Temp1/Temp2)*(.5) 

1940  Temp4>Hatrxd ,1 ,Elllpse*20)*Prob  constant 

1950  Temp5-(Temp4/Temp2)*(.5)  “ 

1 960  Efflin^C  Ellipse )  •-Temp5*Ybar  (Ellipse ) 

1970  Emax_^(Ellipse)-Temp5*Ybar(Ellipse} 

1980  Emin__x(ElliP8e)"-Temp3*Xbar(Ellipse) 

1 990  Emax" x( Ellipse } •Temp3*Xbar ( Ellipse ) 

2000  SUBEND  ~l  End  of  GET  BOUNDS 
2010  I  ~ 

2020  I  CHOOSE  PLOTTER 
2030  1  “ 

2040  SUB  Choose_plotter 

2050  1  Subroutine  to  choose  the  desired  plotter  device 

2060  REPEAT 

2070  Go_on— 1 

2080  DISP  "Do  you  want  to  (D)raw,  (E)rase,  Plot  with  Pen  (0),  (1).  (2)?" 


2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 


LINPUT  Heasagel 
SELECT  Heasagel 
CASE  "0","d" 

PLOTTER  IS  3. "INTERNAL* 
GRAPHICS  ON 
PEN  1 
Go_on«1 
CASE  "O" 

PLOTTER  IS  705,"HPGL" 
GRAPHICS  ON 
PEN  0 
0oj>n-1 
CASE  *1 " 

PLOTTER  IS  705, "HPCL" 
GRAPHICS  ON 
PEN  1 
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2250  COMOn-l 

2260  CASE  “2" 

2270  PLOTTER  IS  705."HPGL" 

2280  GRAPHICS  ON 

2290  PEN  2 

2300  Co  on>l 

2310  CASE 

2320  PLOTTER  IS  3. "INTERNAL" 

2330  GRAPHICS  ON 

2340  PEN  -1 

2350  Go  on-1 

2360  CASE  Use 

2370  PRINT  CHR»(7) 

2380  END  SELECT 

2390  UNTIL  Go  on-1 

21100  SUBEND  ~ 

2410  I 

2420  !  GET  ELL  MEAN 

2430  I  “  “ 

2440  SUB  G«t_«ll_iMan(Elllpae) 

2450  !  This  atTbroutlnoa  proopta  tha  uaer  for  the  mean  of  an  elllpae. 

2460  COM  /Elllpaea/  Matrx(2.2.40).Xbar(20),7har(20),MaJor_axia(20),Minor  ax 
ia(20),Theta(20),Angle(20),Efflln  x(20).Efflax  x(20),Efflin  y(20),Emax  y(20)  ~ 

2470  Clearacreen  “ 

2480  PRINT  TABXX(1.10);"WHAT  IS  X-BAR  for  ELLIPSE  ■; Elllpae s"?"; 

2490  INPUT  Xbar(Elllpae) 

2500  PRINT  Xbar( Elllpae) 

2510  PRINT  TABXY( 1,11) s"HHAT  IS  T-BAR  for  ELLIPSE  ";Elllpaei"7"; 

2520  INPUT  Xbar(Elllpae) 

2530  PRINT  Ybar( Elllpae) 

2540  Pauaeablt 

2550  SUBEND 
2560  1 

2570  1  GET  COVARIANCE 

2580  1 

2590  SUB  Getj:ovarlanoe(Elllpae) 

2600  I  thla  routine  geta  the  covariance  matrix  for  an  elllpae 
2610  COM  /Elllpaea/  Matrx(2,2.40),Xbar(20)tYbar(20) ,Major_axla(20),Mlnor_ax 
la(20),Theta(20),Angle(20),Efflln_x(20),Eoax_x(20),Emln_y(20).Efflax_y(20)  ~ 

2620  Clearacreen  ”  ”  ~ 

2630  FOR  K-1  TO  2 

2640  PRINT  TABXY(1,3*K*6);"ENTER  ELEMENT  (";K;",":K;")  IN  THE  COVARIANCE 

MATRIX  FOR  ELLIPSE*” ;Elllpae;”: 

2650  INPUT  Hatrx(K.K. Elllpae) 

2660  PRINT  Matrx(K.K.Elllpae) 

2670  NEXT  K 

2680  PRINT  TABXY{1,3+9);"ENTER  ELEMENT  (1,2)  IN  THE  COVARIANCE  MATRIX  FOR  E 
LEMENT#";Is": 

2690  INPUT  HatrxC 1,2, Elllpae) 

2700  Hatrx(2,1 , Elllpae )>Matrx(1 ,2, Elllpae) 

2710  PRINT  Hatrxd, 2, Elllpae) 

2720  Pauaeablt 

2730  SUBEND 
2740  I 

2750  I  PAUSEABIT 
2760  I 

2770  SUB  Pauaeablt 

2780  I  pauae  and  wait  for  a  carriage  return 
2790  DISP  "Type  ENTER  to  continue..."; 

2800  INPUT  Garbage! 

2810  SUBBND 
2820  SUB  Clearacreen 
2830  !  Clear  the  acreen 

2840  PRINT  CHR!(12) 

2850  SUBBND 
2860  I 
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2870  i  GET  AXES 
2880  I  “ 

2890  SUB  G«t  ax«s(ElXlpse) 

2900  I  ThTa  subroutine  gets  sn  ellipse  in  terns  of  the  exes  and  the 
2910  I  angle  of  orientation.  These  are  converted  into  a  covariance 

2920  I  natrlx. 

2930  COH  /Ellipses/  Matrx(2,2,ll0).Xbar(20)>ybar(20),Major_axls(20),Mlnor_ax 

ls(20).Theta(20),Angle(20).Ealn_x(20).Emax_x(20),Ealn_jr(20),Eiux_y(20)  ~ 

29*10  Prob_eonstant-2*L0G(.0S)  ~ 

2950  Clearacreen 

2960  PRINT  "UHAT  IS  THE  RADIUS  OF  THE  HAJOR  AXIS  OF  ELLIPSE  "sElllpse;"  :  ” 

s 

2970  INPUT  HaJor_axla( Ellipse) 

2980  PRINT  Major  axla(Elllpse) 

2990  PRINT  *UHAT~ IS  THE  RADIUS  OF  THE  MINOR  AXIS  OF  ELLIPSE  ";Elllpae;''  :  " 

3000  INPUT  Hlnor_axls(Elllpse) 

3010  PRINT  Mlnor~axls( Ellipse) 

3020  PRINT  ”HHAT~ IS  THE  ANCLE  (IN  DECREES)  OF  THE  MAJOR  AXIS  OF  ELLIPSE  ";E 

lllpse:"  :  ”s 

3030  INPUT  Angle(Elllpse) 

30*10  PRINT  AngleC Ellipse) 

3050  i  Convert  to  Covariance  Matrix 

3060  Theta(Elllpse)-PI*Angle(Elllpae)/l80 

3070  Matrxd.l  ,£lllpse)>(MaJor  axis(Elllpse)*C0S(Theta(ElllpBe)))‘‘2*(Mlnor 

axla(Elllpse)<SIN(Theta(Ellipae)))'2  ~ 

3080  Matrx(2,2.£lllpae)>(MaJor  axls(Clllpae)*SIN(Theta(Elllpse))}‘2*(Mlnor 

axia(Elllpse)*C0S(Theta(Elllpse)))*2  ~ 

3090  Hatrx(2,1 ,£lllpae)«((MajQr  axla(Elllpse))*2-(Minor  axls(Elllpse))*2)*C 

OS(Theta(ElHpae))*SIN(Theta(ElllpaeT)  “ 

3100  Matrxd  .2, Elllpae)>Matrx(2,1  .Ellipse) 

3110  FOR  J-1  TO  2 

3120  FOR  K*1  TO  2 

3130  Matrxd,K,Ellip3e)«Matrx(j;K,Elllpse)/(-Prob  consUnt) 

31**0  NEXT  K 

3150  NEXT  J 

3160  PRINT 

3170  Pauseablt 
3180  SUBEND 
3190  ! 

3200  I  DISPLAY  COV 
3210  I  “ 

3220  SUB  01splarjcov(Elllpse} 

3230  I  Displays  the  covariance  natrlx  for  an  ellipse 

32*t0  COM  /Ellipses/  Matrx(2,2,*l0),Xbar(20),Ybar(20),MaJor_axis(20),Mlnor_ax 

ls(20},Theta(20},Angle(20),Enln  x(20),Enax  x(20),Enln  y(20),Enax  y(20) 

3250  FOR  1-1  TO  2  ~  “  ~ 

3260  FOR  J-1  TO  2 

3270  PRINT  Matrxd.J.Elllpse);"  "j 

3280  NEXT  J 

3290  PRINT 

3300  NEXT  I 

3310  SUBEND 

3320  I 

3330  I  ADD  COVARIANCE 
33*10  ! 

3350  SUB  Addj90variance(Flrat, Second, Susner) 

3360  I  This  subroutine  can  be  used  to  add  Sunner-Flrat^'Second 
3370  CON  /Ellipses/  Matrx(2,2.e0),Xbar(20).Ybar(20).Hajor_axla(20).Minor_ax 

la(20).Thata(20).Angle(20).Enln  s(20).£nax  x(20),Eninjr(20).Eaiax  y(20)  ~ 

3380  INTEGER  I,J  “  ~ 

3390  FOR  1-1  TO  2 

3*100  FOR  J-1  TO  2 

3**10  Matrx(I,J,Sunner)-Matrx(I,J.Flrat)*Hatrx(I,J. Second) 

3*120  NEXT  J 

3*130  NEXT  I 
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3«40  SUBEND 
3<I50  I 

3*60  I  COMBINE  ELLIPSE 
3470  I  “ 

34SO  SUB  CoablnejilllpM(Flrst. Second. Conbo) 

3490  !  This  rmclno  finds  the  ^JPL*  eonbination  of  the  First  and  Second 

3500  I  ellipses. 

3510  COM  /Ellipses/  Matrx(2,2.40).Xbar(20).Xbar(20).Najop  axis(20}, Minor  ax 

ls(20),Theu(20).Angle(20),Ealn  x(20).£aax  x(20).EBln  y(20).Egiax  y(20)  ~ 

3520  REAL  Sxl.Sx2  ~ 

3530  I  Calculate  the  new  covariance  natrlx. 

3540  Invert(Flrat,Firat*20) 

3550  Invert(Second,Second*'20) 

3560  Addjeovar lance ( Firat*20 .Second«20 . Conbo*20 ) 

3570  lnv7rt(CoabO'»20,Co(abo) 

3580  I  Calculate  the  new  mean  point. 

3590  Sx1«Matrx(1 ,1  .Flrat«20)*Xbar(Flrat)«Hatrx(1 .2,Flrst*20)*Ibar(Flrst) 

3600  Sxl  aSxl  ^HatrxC  1 .1  ,Second'*'20)*Xbar(Second)4'Hatrx(1  •2,Seeond'»20)*Ibar(Se 

cond) 

3610  Sx2-Matrx(2.1 .  Flrat*'20  )*Xbar(Flrst  )«'Matrx(  2.2,  First«20)*Ybar(  First) 

3620  Sx2>Sx2'»Hatrx(  2 , 1  ,Second'»20 )  *Xbar  ( Second )  ♦Hatrx(  2 , 2 ,  Second*20 )  *Xbar  ( Se 

cond) 

3630  Xbar(Combo)>Matrx(1 .1  .Coabo)*Sx1«'Natrx(1 .2.Conbo)*Sx2 

3640  Xbar  (Combo)  •MatrxC  2,1  ,Cambo)*Sx1«'Hatrx(2,2,Conbo)*Sx2 

3650  SUBEND 
3660  I 

3670  1  OISP  EXTREMES 

3680  I  “ 

3690  SUB  OiapjtxtrenesCElllpsa) 

3700  I  Thl?  routine  displays  the  extreme  x  and  y  values  for  the  specified 
3710  I  ellipse 

3720  COM  /Ellipses/  Matrx(2.2.40),Xbar(20).Ibar(20),Malor  axls(20). Minor  ax 

is(20),Theu(20).An8le(20),Emln  x(20).Eniax  x(20).Emin  y(20),Emax  y(20)  ~ 

3730  PRINT  "E#  ";Elllpaej"T  “ 

3740  PRINT  "Xmin  ":Emin  x(Ellipse); 

3750  PRINT  ”i  Xmax  "tEiax  x( Ellipse); 

3760  PRINT  Ymln  ■  ;Enin3(  Ellipse ) ; 

3770  PRINT  *;  Ymax  ";Emaxy(Elllpse) 

3780  SUBEND 

3790  SUB  DrawjellipseCEllipse, Xmin, naax.Ymin, Ymax, Prob_constant) 

3800  I  ThlT  routine  draws  one  (1)  ellipse  on  the  cuFrent  plotter  device. 
3810  COM  /Ellipses/  Matrx(2,2.40),Xbar(20),Ybar(20).Majop  axls(20) ,Mlnor_ax 

is(20),Theta(20),Angle(20),Enln  x(20).Enax  x(20),Emin  y(20),Eaax  y(20)  ~ 

3820  INTEGER  Sign  ~  “ 

3830  REAL  Xpoint, Upoint.Ypoint, Tempi .Temp2,Temp3,Temp4,Temp5 

3840  GRAPHICS  ON 

3850  I  Draw  top  half  of  the  ellipse 

3860  Signal 

3870  GOSUB  Drawjwlf 

3880  I  Draw  botTom  half  of  the  ellipse 

3890  Sign— 1 

3900  GOSUB  Drawjulf 

3910  I  Finish  uF 

3920  MOVE  (Xbar(Elllpse)-Xmin)*100/(Xmax-Xmin).(Ybar(Ellipse)-Ymin)*100/(Ym 

ax-Ymin) 

3930  DRAM  (Xbar(Ellipse)-Xmin)*100/(Xmax-Xmin).(Ybar(Ellipse)-Ynin)*100/(Ym 

ax-Ymin) 

3940  I  Okey,  dokey,  we're  done  here. 

3950  SUBEXIT 

3960  Draw_half :  I  Draw  half  of  the  ellipse 

3970  If  Signal,  then  draw  top  half;  if  Sign— 1 ,  then  bottom  half. 

3980  FOR  XpolntaEmin  x(Ellipse)  TO  Emax  x(Ellipse)  STEP  .1 

3990  GOSUB  Draw  l^point  “ 

4000  NEXT  Xpoint 

4010  XpolntaEmax  x( Ell ipse) 

4020  GOSUB  Draw  T  point 
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0030  RETURN 
«0«0  I 

4050  I  Coopute  ecch  point  «nd  draw  the  new  line.  (It's  here  since  we  cell 

4060  I  it  twice 

4070  1 

4080  Drsw_1_peint:  I 

4090  '*Upoint>Xpoint-Xbsr(Ellipae) 

4100  Temp1>Hstrx(1 ,2.Elllpae'»20)*Upoint 

41 1 0  TeMp2-Temp1 *Temp1 -Hatrx(2 .2.Ellipse*20 )*(Mstrx(  1.1. Ellip3e*20 )*Upoint ■ 

Upoint^Prob  constant) 

4120  I?  Temp2<10‘(-10)  THEN  TeBp2-0 

41 30  Xpoint-< -Teopl *Sign*SQR(Temp2 ) ) /HatrxC  2 .2 , Elllpse«20 ) *Ybar (Ellipse ) 

4140  IF  Xpoint<-Efflin  x( Ellipse)*. 001  THEN 

4150  MOVE  (Xpoint~Xnin)*100/(XMX-Xmin),(Ypoint-Xain)*100/(Yaax-Xnin) 

4160  ELSE 

4170  DRAM  (Xpoint>Xmin)>t00/(XBax-XBin).(Ypoint-Yain)*100/(Xaax-Tmin) 

4160  END  IF 

4190  RETURN 

4200  SUBEMD 

4210  I 

4220  I  TEST 

4230  ! 

4240  SUB  TestCFirst, Second, Work. Teat  val) 

4250  I  nils  routines  calculates  the  acceptance  test  criteria  for  First  and 

4260  I  Second.  Work  la  used  as  a  "scratchpad”. 

4270  COM  /Ellipses/  Matrx(2.2.40).Xbar(20).Ybar(20).MaJor_axis(20),Minor_ax 

ls(20),Theta(20),Angla(20).Eain_x(20),Eoax  x(20).EBln_jr(20),EBSx_y(20)  ~ 

4280  Addjsovar lance (First .Second . HorkT 

4290  Inv?rt(Work,Hork*20) 

4300  0iffx>Xbar(Firat)-Xbar(Seeond) 

4310  Dlffy-ybar(Firat)-Ybar(Seeond) 

4320  Teapl-Matrxd ,1 ,Hork»20)»Dlffx*Matrx(1 ,2,Hork*20)»Diffy 

4330  Teap2-Matrx(2,1 ,Hork*20)«Dlffx*Matrx(2,2,Hork»20)»Diffy 

4340  Teat  val-01ffx«Tei«pl*Dlffy*Temp2 

4350  SUBEND  " 

4360  SUB  Sinjtllipse 

4370  1  Ellipse  Conbinatlon  Sisulation  Program 

4380  1  1/16/85  Update  7/18/85 

4390  ! 

4450  I  Specify  the  cooaon  variables 

4460  I 

4470  COM  /Ellipses/  Matrx(2,2,40).Xbar(20),Ybar(20),Major_axls(20),Minor_axi 

a(  20 )  .Theta  ( 20 ) ,  Angle(20 ) ,  Emin_x(  20 ) ,  Eaax_x(  20 ) ,  EBin_y  ( 20 ) ,  Eaax_jr(  20 ) 

4480  COM  /Driver/  0ba(2)  ~  “ 

4490  I  Variable  Definitions 

4500  I 

4510  INTEGER  EllnuB,Iaat.XBin,Xaax.Yain.YBax 

4520  INTEGER  Nobaervatlona(2) 

4530  DIM  Sx(2,10),Xa(5000).Ts(5000) 

4560  I 

4570  I  Use  of  HATRX  Array 

4580  I 

4590  I  MATRIX  1  is  an  initial  aatrlx 

4600  I  MATRIX  2  is  an  initial  Matrix 

4610  I  MATRIX  3  is  the  theoretical  coablnation  of  8  and  9 

4620  I  MATRIX  4  la  the  eatiaate  of  8 

4630  I  MATRIX  5  is  the  est lasts  of  9 

4640  I  MATRIX  6  is  the  coablnation  of  4  and  5 

4650  I  MATRIX  7  is  the  aua  of  3  and  4  (for  test) 

4652  I  Note  that  if  the  True  Variance-Covariance  is  used  instead  of  the 

4653  I  the  eatiaate.  this  is  equal  to  the  sua  of  8  and  9. 

4654  1  MATRIX  8  is  MATRIX  1  divided  by  sample  size 

4655  I  MATRIX  9  is  MATRIX  2  divided  by  sample  size 

4660  I  MATRIX  21-29  are  the  inverses  of  1-9 

4670  I 

4680  I  Initialization 
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11690  I 

4700  GOSUB  Init  aubi 

*710  !  ~ 

*720  I  Main  Prograa  Loop 

*730  I 

*7*0  Naln_loop:  I 

*750  ~ClaaracrMn 

*760  raiMT  "Elllpaa  SlBulatlona* 

*770  PRINT 

*780  PRINT  "Entor  the  Latcar  of  your  cholca:" 

*790  PRINT  "  (A)  Raaat  tha  Random  Number  Seed* 

*800  PRINT  *  (B)  Simulate  uaing  true  Covariance  Matrleea" 

*810  PRINT  *  (C)  Simulate  uains  estimated  Covariance  Matrices" 

*820  PRINT  ■  (0)  Display  Results  on  the  Printer" 

*830  PRINT  "  (X)  Exit  Prograa" 

*831  Get_optlon!  ! 

*8*0  INPUT  "Enter  Your  Choleet*.Malncholee$ 

*860  SELECT  MalneholceS 

*870  CASE  "A", "a" 

*880  GOSUB  Get  seed 

*890  CASE  ■B","b^ 

*900  Use  trua>l 

*910  GO^B  Simulate 

*920  CASE  «C","c" 

*930  Use  true-0 

*9*0  GO^B  Simulate 

*950  CASE  "0",*d" 

*960  GOSUB  Display  results 

*970  CASE  "X","x"  " 

*980  GOTO  End  subi 

*990  CASE  ELSE  ” 

5000  PRINT  CHR*(7) 

5010  GOTO  Get  option 

5020  END  SELECT 

5030  GOTO  Main  loop 

50*0  ! - - 

5050  I  Utility  Subroutines 

5060  1 - 

5070  1 

5080  I  Initialize  Prograa 

5090  I 

5100  Inlt  subi:  ! 

5101  “randomize 

5110  Prob  constanta2*L0G(.05}  I  Confidence  level  parameter 

5111  FOR  T-1  TO  2 

5112  FOR  J-1  TO  2 

5113  Matrx(I,J.8)-Matrx(l.j.i)/0bs(1) 

51 1  *  Matrx( I , J . 9 )-Hatrx( I , J .2 ) /Obs( 2 ) 

5115  NEXT  J 

5116  NEXT  I 

5117  Xbar(8)-Xbar(1) 

5118  Ibar(8)-tbar(1) 

5119  Xbar(9)-Xbar(2) 

5120  Ybar(9)-Ybar(2) 

5122  I  Set  parameters 

5230  RETURN 

53*0  I 

5350  I  Generste  the  new  ellipses  using  random  observations 

5360  I 

5370  Generate  new:  l•••PLAG••• 

5380  FOR  Iaat-1  TO  2 

5390  M-Oba(Iaat) 

5*00  I  Generate  the  X  and  Y  values 

5*10  Xsua^ 

5*20  Ysum-O 

5*30  Ce>SQR(Matrx(l,l.Imat)) 


5<i40  A-Slgu(1.2.lMt)/Cc 

5450  Be€-SQR(Hatrx(2,2.Imat)-A*A) 

5460  FOR  1-1  TO  M 

5470  U1-RND 

5480  U2-Rin) 

5490  X1-SQR{-2»L0C(01 ))»SIM{2»PI»U2) 

5500  Y1 -SQR( “2*L0C{ U2 ) ) •C0S(2«PI«U1 3 

5510  Xs(I)-Cc«X1-Xb«r(lBat) 

5520  Ts(I)-A»X1-BM*yi*Xbar(Ifflat) 

5530  Xsua-3Caua«Xa  ( 1 3 

5540  Ysun-Xaua-XaCli 

5550  NEXT  I 

5560  3Cbar(Imat-3)-Xauai/H 

5570  Xbar(laat*3)-TauB/M 

5580  I  Calculate  the  varlaneea 

5590  Hatred  ,1  (laat-3)~0 

5600  Hatrx(2,2,Iaat-33-0 

5610  Hatred  .2. laat-33-0 

5620  FOR  1-1  TO  H 

5630  Hatrxd.1.lBat*33-Hatrxd.1,lBat*33*(Xa(l3-Xbar(lBat-3)>*(Xa(l3-Xba 

r(Iaat*3)3 

5640  Hatrx(2,2,lBat*33-Hatrx(2.2.I«at*33»(Ya(I3-Ybar(I*at*3)3»(ya{I3-yba 

r(lBac-333 

5650  Hatrxd.2,lBat-3)-Hatrxd.2,lBat*33*(X8(I3-Xbar(Iinat*333»(Ya(I3-Yba 

r(Iaat*333 

5660  NEXT  I 

5670  Hatred  .1  .lBat*33-Hatred  .1  ,lBat*33/((H-l  )»M3 

5680  Hatrx(2,2.lBat'»33-Hatrx(2,2.lBat*3>/((H>1 3*M3 

5690  Hatred  .2, :Bat-33-Hatred  ,2,lBat*33/((H-1 3»M3 

5700  Hatrx(2,1  ,XBat*3)-Hatrxd  .2.1aat'»3} 

571 0  NEXT  taat 

5720  RETURN  !  Generate  new 

5730  t 

5740  t  Generate  the  inveraea  of  our  matricea 

5750  I 

5760  Gen  Inveraea i  ! 

5830  Patchit  t 

5840  IF  Uae_true-0  THEN  Patch2 

5850  I  To  uae  the  true  oatrlcea,  copy  thea  fron  8  and  9  Into  4  and  5 

5860  FOR  lut-8  TO  9 

5870  FOR  1-1  TO  2 

5880  FOR  J-1  TO  2 

5890  Hatrx( i . j . lmat-4  3-Hatrx( I . J . loat  3 

5900  NEXT  J 

5910  NEXT  I 

5920  NEXT  lut 

5930  Patch2:  I 

6151  CoBblnejilllpae(8,9.3) 

6152  Coablne'~elllpae<4,5.6) 

6170  RETURN  I  Gen  Inveraea 

6430  I - = - 1 

6440  I  Prograa  Subroutlnea  I 

6450  1 - ! 

6460  I 

6470  I  GET  SEED  (Henu  Henu  Choice  ’•C''3 

6480  I 

6490  Get_aeed :  I 

6500  ~Cle8r8creen 

6510  PRINT  "ENTER  A  SEED  (1  to  2*31-23:  "j 

6520  INPUT  Seed 

6530  Seed4-VAL$(Seed3 

6540  PRINT  Seeds 

6550  RANDOHIZE  Seed 

6560  PRINT 

6570  Pauaeabit 

6580  RETURN 
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6590  I 

6600  i  SIMUUTE 

6610  I 

6620  SiBulate :  I 


6630 

6640 

UH? 

6650 

6660 

6670 

6680 

6690 

6700 

6710 

6720 

6730 

6740 

6750 

6760 

6770 

6780 

6790 

6800 

6810 

6820 

6830 

6840 

6860 

6860 

6870 

6880 

6890 

6900 

6910 

6940 

8))'2 

6950 

6960 

9))‘2 

6970 

6971 
6980 
6990 
7000 
7010 
7020 
7030 
7040 
7050 
7060 
7070 
7080 
7090 
7100 
7110 
7120 
7130 
71 40 
7150 
7160 
7170 
7180 
7190 
7200 
7210 
7220 


Cl««rserMn 
PRINT  TABXXd.lT)!’ 


HOU  HANT  SIMULATIONS  DO  XOU  WISH  TO  R 


INPUT  Nalaulatlona 

PRINT  Nalfflulatlona 

H_accepe«0 

Njrajact«0 

Anotlnl-O 

Anotln2aO 

Anotlnb^O 

Rnotlnl >0 

Rnotln2-0 

Rnotlnb-O 

Pauaaablt 

Claaraoraan 

PRINT  TABXT(5.3)t*NUHBCR  OF  SIMULATIONS*  0" 

PRINT  taB«(5, 4) 

PRINT  TABXY(5.5)s”ACCEPTE0  ERROR  ELLIPSES-  0* 

PRINT  TABXT(5.6)t”REJECTEO  ERROR  ELLIPSES-  0- 
PRINT  TABXY(5.9)j"ACCEPTED  BREAKOUT" 

PRINT  TABXY(5.10);"Targat  Qn«  la  NOT  In  Coabinad  Ellipae-  0" 

PRINT  TABXT(5.11):"Target  Two  la  NOT  in  Coabinad  Elllpaa-  0" 

PRINT  TABXT(5.12):"Nalthar  Targat  la  In  Coabinad  Elllpaa-  0” 

PRINT  TABXX(5,1*):”REJECTE0  BREAKOUT" 

PRINT  TABXY(5.15);"Targat  Ona  la  NOT  in  Coabinad  Elllpaa-  0" 

PRINT  TABXY(5.l6);"Targat  Two  la  NOT  in  Coabinad  Elllpaa-  0" 

PRINT  TABXT(5,17);"R*lther  Targat  la  In  Coabinad  Elllpaa-  0" 

FOR  Iala-1  TO  Nalaulatlona 
GOSUB  Ganarata_naw 
GOSUB  Gan  invaraaa 

Ona_in-Matrx(  1 , 1 , 26 )  •  (Xbar  ( 6 ) -Xbar  ( 8 ) )  *'2*Hatrx  ( 2 , 2 , 26 )  •  ( Tbar  ( 6 ) -Ybar  ( 

Ona  In-Ona  ln*2»Matrx<1 ,2,26)«(Xbar(6)-Xbar(8))»(Ybar(6)-Ybap(8)) 
Two“ln-Matrx(1.1,26)»<Xbap(6)-Xbar(9))‘2»Matrx(2,2,26)*(Ybap(6)-Ybar( 

Two  In-Two  ln*2*Matrx(1 ,2,26)»(Xbar(6)-Xbar(9))»(Ybar(6)-Ybar(9)) 
Taa7(4,5,77Taat2) 

! 

1 

PRINT  TABXY(28.3);Iala 
IF  Ta8t2<-Prob_conatant  THEN 
I  Aecapt  aa  Maa 
N  aecapt-N  accapt*l 
PlflNT  TABxT(30.5)jN_accapt 
IF  Ona_in>-ProbjBonatant  THEN 
AnotTnl  -Anot  lITl  - 1 
PRINT  TABXY(44,10)tAnotln1 
END  IF 

IF  Two_ln>-Probjeonatant  THEN 
AnotTn2 - Ano t ln2 * 1 
PRINT  TABXY(44.11):Anotin2 
END  IF 

IF  Ona_In>-ProbjBonatant  AND  TVo_ln>-Prob_c«matant  THEN 
AnocTnb-Anot inb  * i 
PRINT  TABXY(44.12)tAnatlnb 
END  IF 
ELSE 

N  raJaet-N  rajact-t 
PTTINT  TABXT(30,6)!N_pajact 
IF  Ona  in>-Probjeonatant  THEN 
RnotTnl  -Rnotlril  -1 
PRINT  TABXY(44,15);Rnotlnl 
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7230  END  IF 

7240  IF  Two  ln>-Frob  constant  THEN 

7250  RnotTn2>linotin2'»l 

7260  print  TABXX(44,l6)jRnotln2 

7270  END  IF 

7280  IF  Ons_tn>-Prob_eonatant  AND  Two  ln>-Prob  constant  THEN 

7290  RnotlnbaRnotinb*!  “  ~ 

7300  print  TABXT(44,17)sRnotlnb 

7310  END  IF 

7320  END  IF 

7330  NEXT  Islm 

7340  PRINT  TABXY(5.19);"»»*  SIMULATION  COMPLETE  •••"j 

7341  BEEP 

7342  BEEP 

73«3  BEEP 

7344  BEEP 

7350  Pausoabit 

7360  RETURN 

7370  ! 

7380  I  Display  results 

7390  I 

7400  Dlsplay_results : ! 

7410  Clearscreen 

7420  PRINTER  IS  9 i WIDTH  132 

7421  PRINT 

7422  PRINT 

7423  PRINT 

7424  PRINT 

7430  PRINT  "Results: 

7440  IF  Use_true-1  THEN 

7450  PRINT  "Using  TRUE  Covariance  Matrices" 

7460  ELSE 

7470  PRINT  "Using  ESTIMATED  Covariance  Matrices" 

7480  END  IF 

7490  PRINT 

7491  FOR  loat-l  TO  2 

7492  PRINT  "Base  Distribution  #";Iaat 

7493  PRINT  "Mean:  (":Xbar(Isiat):",  ";Tbar(lBat) 

7494  PRINT  "Observations:  "sObsdmat) 

7495  PRINT  "Covariance  Matrix:" 

7496  PRINT  "*  ";Matrx(1 ,1  .Isiat)  s"  ";Matrx(1 ,2,Iiiist) ;"  *" 

7497  PRINT  "•  ";Matrx(2,1  ,Imat) ;"  "sMatrx(2,2,Iniat)  j"  •" 

7498  PRINT 

7499  NEXT  Inst 

7500  FOR  lMf4  TO  5 

7501  PRINT  "Ellipse  #";I«Mit-3 

7502  PRINT  "Mean:  (":Xbar(Iaat) "sTbardmat) 

7504  PRINT  "Covariance  Matrix:" 

7505  PRINT  "•  "sMatrxi  1.1, last)  t"  ";Matrx(1 ,2.Iaiat) ;"  •" 

7506  PRINT  "•  "5Matrx(2,1,Iaat){"  ";Matrx(2,2,lBat) ;"  •" 

7507  PRINT 

7508  NEXT  iMt 

7590  PRINT  "Last  Coablned  Ellipae:" 

7600  PRINT  "Moan:  ("jXbar(6)j",  "iXbar(6)i") ;" 

7610  PRINT  "Covarlanoe  Matrix:" 

7620  PRINT  "•  "iMatrxd,  1.6):"  ";Matrx(1 ,2,6) ;"  •" 

7630  PRINT  "*  ■|Natrx(2,1,6)s"  ";Matrx(2,2,6) :"  •" 

7640  PRINT 

7650  PRINT  "Slaulatlon  Results:" 

7660  PRINT  "«  of  slsulations:  ";NalMlations 

7661  PRINT 

7662  PRINT  "  Totals:  %  of  Total:  %  of  Category" 

7670  PRINT  "Accepted:  ":N  accept:"  ":100"N  accept/Nainulationa 

7680  PRINT  "  Target  1  NOT  in:  ”:Anotlni:  “ 

7681  PRINT  "  ":100*Anotin1/NaiMlations:"  ":l00*Anotin1/N  accept 

7690  PRINT  "  Target  2  NOT  in:  "lAnotinE;  “ 


7691  PRIMT  "  ■;100»Anotln2/Haimulatlona5" 

7700  PRINT  "  Naithar  ona  In:  ";Anotlnbi 

7701  PRINT  »  ";100«Anotlnb/Nai«»latlona;" 

7710  PRINT  "Rajactad:  •;N_paject;" 

7720  PRINT  "  Targat  1  NOT  in:  "iRnotlnl; 

7721  PRINT  "  ":100*Rnotlnt/Naiaiulatlona;* 

7730  PRINT  "  Targat  2  NOT  In:  •;Hnotln2; 

7731  PRINT  "  •;lOO*Rnotin2/Naimulationa:” 

77*10  PRINT  "  Naithar  one  in:  *;Rnotinb 

77*11  PRINT  "  ";100*Rnotlnb/NaiBUlationa:'' 

7750  PRINT 

7751  PRINT 

7753  PRINT  CHR$(12) 

7760  PRINTER  IS  1 t WIDTH  80 

7770  RETURN 
7780  End  aubi:! 

7790  “SUBIND 

7800  SUB  PowarjaUlpaa 

7810  I  Targat  Ullpaa  Chl-aquara  Taat  Powar  Calculation 
7820  I  Specify  coonon  variablea 
7830  COH  /Elllpaea/  Hatrx(2.2.*i0).Xbar(20),1 


■ ;  1 00  * Anot In2/Njtecep  t 

■ ;  100*Anotlnb/N_accept 
* ;  1 00*N_raJect/Naionilatlona 

•j i00»Rnotln1 /N_reject 

■ j 1 00»Rnotln2/N_reJect 

" ;  100**  Rnotinb/N_r  eject 


7bar(20) ,HaJor_azia(20) ,Minor_axi 
a(20),Theta(20).Angle(20).Eiiiin_x(20),Emax_x(20).Efflln jr(20),Eiiiax_y(20) 


78*10  COH  /Driver/  0ba(2) 

7850  I  Variable  Definition 

7870  DIM  Newjnat(2.2.2).Delta(2) 

7880  DIM  Table(62.2) 

7900  !  Initialize  Program 

7910  GOSUB  Inlt  aub2 

7920  I  MENU  LOOP 


7930  MenuJLoop: 


1 


79*10  Clear_acreen 
7950  GOSUB” Power  atuff 

7960  SUBEXIT  “ 

7970  1 

7980  t  INITIALIZE  THE  PROGRAM 

7990  I 

8000  Init_aub2:  I 

8010  Prob_conatant-2*L0G( .05 ) 

8020  I  Load  in  the  Power  Table 

8030  RESTORE  Powar  data 

80*0  FOR  1-1  TO  62~ 

8050  FOR  J-1  TO  2 

8060  READ  Tabla(I.J) 

8070  NEXT  J 

8080  NEXT  I 

8170  RETURN 
8180  I 

8190  I  Powar  calculationa:  find  the  non-centrality  paraMtar,  Lambda. 

8200  I  to  uaa  with  the  tablaa. 

8210  ! 


8220  Powarjituf f : I 

8230  p7inT  "Power  Calculationa:  Non-centrality  Paraaatar" 

8240  PRINT 

8250  FOR  1-1  TO  2 

8260  FOR  J-1  TO  2 

8270  Powar  matCI.J.I )-Hatrx(I,j,1 )/0ba(l )*Matrx( I.J.2)/0ba(2) 

8280  NEXT  J  ~ 

8290  NEXT  I 

8300  I compute  inveraa  of  powar  jaat 

8310  Det-PowarjaatCI  ,1,1  )**PowaF  mat(2,2,1  )-Powar_mat(1 ,2,1  )"Powar_mat(1 ,2,1 ) 
8320  Powar_matTl  ,1 ,2)-Powar_matT2,2.1 )/Dat  ” 

8330  Power3d4(2,2,2)-PowarjBat(1 ,1 ,1  )/Dat 

83*0  Powarjaatd  ,2,2)— Powar_mat{1 ,2.1  )/Det 

8350  Powar'nat(2,1 ,2)— Powar“mat(2,l  ,1  >/Dat 

8360  print" "Enter  the  diffarFneaa  in  the  Mean  eomponanta:" 

8370  PRINT 
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6380  PRINT  "Mean  difference  along  the  X  axis: 

8390  INPUT  "X  Difference :",Delta(1 ) 

81(00  PRINT  DelUd) 

8(I10  PRINT  "Mean  difference  along  the  X  axis:  "; 

81(20  INPUT  "Y  Difference:". 0elta<2) 

8((30  PRINT  Delta(2) 

8(|(I0  i  Calculate  parasMters 

8((50  Lambda-0 

8460  FOR  1-1  TO  2 

8470  FOR  J-1  TO  2 

8460  Lambda-Lambda*Powerjiiat(I,J,2)*Delta(I)*Delta(J) 

8490  NEXT  J 

8500  NEXT  I 

8510  PRINT  "Lambda  -  "iLambda 

8520  PRINT 

8530  I  Get  Linear  Interpolation  Result 

8540  GOSUB  Linear_lnterp 

8550  PRINT  "Linear  Power:  ";Lijpower 

8560  I  Get  Lagrange  Interpolation  Result 

8570  GOSUB  Lagrange_lnterp 

8580  PRINT  "Lagrange  Power:  ":La_power 

8590  Pauseablt 

8600  IF  Delta(1)<>0  OR  0elta(2)<>0  THEN  8360 

8610  RETURN 
8620  I 

8630  I  Linear  Interpolation 

8640  I 

8650  Linear_interp :  I 

8660  I  Find  Bounding  Values 

8670  Search-I 

8680  IF  Umbda>-39  THEN 

8690  Li  power-1 

8700  return 

8710  END  IF 

8720  WHILE  Umbda>Table(Search.1 ) 

6730  Search-Search* 1 

8740  END  WHILE 

8750  IF  Lambda-Table(Search,1 )  THEN 

8760  Li  power-Table(Search,2) 

8770  R^RN 

8780  END  IF 

8790  Lambda_h-Table(Search,1 ) 

8600  Lambda~l*1'Abl*(S«4''ch~'' d ) 

6810  Power_jr-Table(Search,2) 

6620  Power^'fableCSearch-l  ,2) 

8830  IF  Power_h-Power_l  THEN 

8840  Li  power-Power~h 

8850  return  ~ 

8860  END  IF 

8870  Powerl  ■Powor_h*(Lambda-LagAda_l)/(Lambda_h-Lambda_l) 

8880  Power2-Power~l*  ( Lambda_h-Lambda)  /  ( Lambda~h-Lambda~l ) 

8890  Ll_power-Power1*Power2~ 

8900  RETURN 

8910  I  Lagrange  Interpolation 

8920  Lagrange_interp :  I 

8930  Sear?h-i 

8940  IF  Laabda>-39  THEN 

6950  La  power-1 

8960  return 

8970  END  IF 

8960  WHILE  Unbda>Table(Search.1) 

8990  Search-Search* 1 

9000  END  WHILE 

9010  IF  Unbda-Table(Search,l)  THEN 

9020  La  power-Table(Search,2) 

9030  RWRN 
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9050  Suoaer'O 

9060  FOR  I-Seareh-3  TO  Searcn*2 

9070  Prod«Table(I,2) 

9080  FOR  J-Search-3  TO  Search*2 

9090  IF  JOI  THEN 

9100  Prod-Pro<l«(Lambda-Table(J.1))/(Table(I,1  )-Table(J,1 )) 

9110  END  IF 

9120  NEXT  J 

9130  Suomer-SuiiiDer*Prod 

91 «0  NEXT  I 

91 80  Ua_power»Sufflaer 

9160  RETURN 

9170  IChi-9quare(2)  Power  Table,  lat  coluon  la  non-centrality  parameter, 

9180  12nd  la  power.  0.05  algnlflcance  label.  From  Selected  Tablea  in 

9190  IMathematlcal  Statlatica,  Voluaie  1. 

9200  Power_data :  I 
9210  DATA  .0..05 

9220  DATA  .1..0576 

9230  DATA  .2,. 0653 

92*10  DATA  .3.. 0733 

9250  DATA  .<(..0810 

9260  DATA  .5.. 0896 

9270  DATA  .6,. 0980 

9280  DATA  .7,. 1065 

9290  DATA  .8. .1151 

9300  DATA  .9. .1239 

9310  DATA  t..1327 

9320  DATA  1.2,. 1507 

9330  DATA  1.9, .1691 

93“0  DATA  1.6,. 1877 

9350  DATA  1.8,. 2065 

9360  DATA  2.0,. 2255 

9370  DATA  2. 2,. 2997 

9380  DATA  2. 9,. 2639 

9390  DATA  2. 6,. 2831 

9900  DATA  2. 8,. 3029 

9910  DATA  3.0, .3215 

9920  DATA  3. 5,. 3690 

9930  DATA  9,. 9159 

9990  DATA  9. 5,. 9609 

9950  DATA  5,. 5037 

9960  DATA  6,. 5890 

9970  DATA  7,. 6559 

9980  DATA  8,. 71 76 

9990  DATA  9,. 7707 

9500  DATA  10,. 81 59 

9510  DATA  11,. 8526 

9520  DATA  12,. 8832 

9530  DATA  13,-9080 

9590  DATA  19,-9280 

9550  DATA  15,-9990 

9560  DATA  16,-9567 

9570  DATA  17.. 9667 

9580  DATA  18,-9795 

9590  DATA  19,-9805 

9600  DATA  20,. 9852 

9610  DATA  21.. 9888 

9620  DATA  22.. 99 16 

9630  DATA  23,-9937 

9690  DATA  29,-9953 

9650  DATA  25.. 9965 

9660  DATA  26,-9979 

9670  DATA  27,-9981 

9680  DATA  28,-9986 

9690  DATA  29.-9989 
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VV 


9700 

DATA 

30.. 9992 

9710 

DATA 

31..  9999 

9720 

DATA 

32.. 9996 

9730 

DATA 

33.-9997 

97«0 

DATA 

39.. 9998 

9750 

DATA 

35.. 9998 

9760 

DATA 

36.-9999 

9770 

DATA 

37.-9999 

9780 

DATA 

38.. 9999 

9790 

DATA 

39.1.0000 

9800 

DATA 

90.1.0000 

9810 

DATA 

91 . 1 .0000 

9820 

DATA 

92.1.0000 

9830 

End  sub2 :  1 

98X0 

SUBEMD 

C.  ELLIPSTUFF;  Ellipse  Routine  Library 


Elllpstuff  is  an  "ellipse  routine  library,"  vrtiich  contains  a  large  number  of 
routines  to  make  working  with  ellipses  easier.  As  such,  both  Genellipse  and 
Ellsim  contain  all  of  these  routines,  but  they  are  also  listed  separately. 

If  more  programs  need  to  be  written,  these  are  the  routines  to  build  them 
with. 


i 


10 

11 

12 

13 

20 

30 

31 

32 
HO 
50 
60 
70 
80 
81 
90 
100 
110 
1000 
1010 
1020 
1030 
lODO 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 

lino 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

i2no 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

13«0 

1350 

1360 

1370 

1380 

1390 

moo 

mio 

1420 

1*30 

1*31 

1*32 

1*33 

1*3* 


I  TEST  PROGRAM  FOR  ELLIPSE  ROUTINE  LIBRARY 
I  LIBRARY  SUBROUTINES  BEGIN  ON  LINE  1000. 

t  THIS  TEST  PROGRAM  GETS  AN  ELLIPSE  FROM  THE  USER  AND  OISPUYS  IT  ON 
I  THE  SCREEN. 

GCLEAR 
GRAPHICS  ON 
REAL  ProbjeonstMt 
Prob_eona?anca2*L0G( .05 ) 

G«tjill_0M«n(1 ) 

I  C«t~covirlanc«(1 ) 

Cbt^axMd  } 

Inv7rt(1 .21 ) 

CatJbounda( 1 .Probjeonstant) 

Oidose _plott»r  ~ 

Draw  alllpaad  .■>50.50,-50.50.Prob  constant) 
stop"  " 

END 

t  ELLIPSE  ROUTINE  LIBRARY.... 

1  This  flic  contains  standard  alllpso  routines  and  data  variables, 

I  including  tbe  following: 

I  —  Entry  of  ellipses  by  covariance  Batrices. 

I  —  Entry  of  ellipses  by  axes  and  orientation. 

I  Display  of  ellipses  on  screen  and  plotter. 

I  2x2  Matrix  inversion  routine  for  use  with  Matrx. 

I  —  Matrx,  an  array  which  stores  covariance  aatrices  and  their 
I  inverses  (up  to  10  natrlces). 

1  —  Axes/orientation  to  Covariance  natrlx  conversion  routine. 

! 

1  USING  MATRX 

1  Matrx  is  designed  to  hold  covariance  matrices  and  their  inverses 
!  for  10  ellipses.  In  general,  Matrxd )...Matrx(20)  are  the  covariance 
!  matrices  and  Matrx(2l )...Matrx(*0)  are  the  corresponding  Inverses. 

I  This  is  the  convention  assumed  by  a  number  of  these  routines. 

I  In  oases  where  the  inverse  is  calculate  first,  put  in  Matrx(25),  say, 

I  and  then  call  Invert(2S.  5).  This  will  put  the  covariance  matrix  in 
t  Matrx(5). 

I 

I  THE  ROUTINES  ARE  CALLED  AS  FOLLOWS: 

I  Invert(Srca,Oest) 

I  This  oomsund  will  invert  the  covariance  matrix  in  Matrx(Sren)  and 
I  put  the  result  in  Hatrx(Dest) 
i  GetJBoundsCEllipse.Probjeonstant} 

I  This  command  will  geT  the  extreme  points  of  the  ellipse.  Note  that 
I  the  ellipse's  covariance  matrix  must  have  been  inverted. 

I  Draw^Ellipae ( Ell ipse , Xmin , xmax , Ymin , Xmas , Probjsonstant ) 

I  This  will  draw  the  given  ellipse.  The  remaining  variables  define 

I  the  screen.  Note  that  Get  bounds  must  have  been  executed. 

I  Get_EllJ1ean(  Ellipse)  " 

I  %is" is  an  input  routine  to  read  in  Xbar  and  Ybar  for  the  given 

I  ellipse. 

I  Get  CovarianoeCEllipse) 

I  This  is  an  input  routine  to  read  in  the  covariance  matrix  for  the 
I  given  ellipse. 

I  Get_Ases( Ell ipse) 

I  This  is  an  input  routine  which  reads  in  the  axes  and  orientation 
I  of  the  ellipse,  and  converts  this  to  covariance  matrix  form. 

I  Choose_Plottor 

I  Thi7  prompts  the  user  to  draw  the  ellipse  on  the  plotter  or  the 

I  screen. 

1  Dlaplay_Cov(Ellipse) 

I  niis"oommand  displays  the  given  covariance  matrix. 

I  Teat(Ellipse1 ,Ellipae2.Hork.TestJStat) 

I  This  caloulatss  the  acceptance  teat  criteria  for  the  two  ellipses. 

I  Work  is  a  Matrx  entry  used  as  working  apace. 

I  Add  Covarianee(ElliDsel.Elllose2. Sumer) 
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•  •*■*4  •*  4** 

<435  i  This  adds  any  tuo  Natrx  antrlea  Into  a  third  Katrx  entry. 

1436  i  CoablneJElllpaa(Elllpae1.EUlpse2.Conbo) 

1437  I  This  coablnea  any  two  Hatrs  entires  Into  a  third  Matrx  entry. 

1438  I  01sp_Extreawa(Elllpse) 

1439  I  TiTls  displays  the  extreae  X  and  T  bounds  for  the  ellipse. 

1441  I  Also  Included  are  Pauaeablt  and  Clearsereen. 

1450  1 

1 460  I I 
1470  SUB  InvertCSrce.Dest) 

1480  I  This  routine  inverts  any  covariance  aatrlx  In  Hatrx  and  places 
1490  I  the  Inverted  matrix  In  Beat. 

1500  COM  /Ellipses/  Katrx(2.2,40),Xbar(20),Xbar(20).Major_axls(20),Mlnor_ax 

la(20),Thata(20).Angle<20).Ealn_x(20),EBax_x(20).Enlnjr(20),Enax_)r(20)  ~ 

1510  Det-Hatrx(1 ,1  .Srce)*ltatrs(2.2,SrM>-Matrx(1 .2,Sree)*Matrx(1 .2.Srce) 

1520  Hatrxd  ,1 ,0eat)«Hatrx(2.2,Srce)/Det 

1530  Hatrx(2,2.0est)aMatrx(1 .1 .Srce)/Dst 

1540  Hatrxd  ,2,Desc)>-Hatrxd,2,Sree)/Det 

1550  Hatrx(2.l ,0eat)>-Matrx(2,1.Sree)/Det 

1560  SUBENO  I  End  of  SUB  Invert 

1600  SUB  Get  boundsCElllpse.Prob  constant) 

1610  I  Tills  subroutine  ealcuTates  the  X  and  T  limits  for  the  given  ellipse 
« 

1620  COH  /Ellipses/  Matrx(2,2.40),Xbar(20),Xbar(20),Hajor_axis(20).Mlnor_ax 

la(20),Theta(20).Angle(20),Ealn_x(20),£aax_x(20),Enin _p(20),Eiux _jr(20)  ~ 

1630  REAL  Tempi ,TeHp2.TenpTtTemp4, Temps 

1 650  Tempi  •Hatrx(  2.2.  Elllpae'»20 )  *Prob_constant 

1660  Temp2-Hatrxd  .2.Elllpse«20)*Matrx(1 .2,Elllpae*20)-Matrxd  .1  .Elllpse*20 

) *Ma trx( 2 . 2 . Ell 1 pae 420 ) 

1670  Temp3*(Temp1/Temp2)*(.5) 

1680  Temp4-Hatrxd  ,1  ,Elllpae420)*Probjeonstant 

1690  Temp54(Temp4/Tamp2]*(.5) 

1 700  Emlnjf  ( Ellipse  )<i-Teap54Tbar<Eillpse ) 

1710  Emax^(Elllpse)«TempS4Xbar(Elllpse) 

1 720  Emln_x(  Ellipse  )»Teap3‘»Xbsr(  Ellipse ) 

1730  Emax~x(Elllpae)-Temp3*Xbar(Elllpse) 

1740  SUBEND  "l  End  of  GET_B0UNDS 
1780  SUB  Choose_plott4r  ** 

1790  I  Subroutine  to  choose  the  desired  plotter  device  (also  contained 

1791  I  in  GRAPHER). 

1800  REPEAT 

1810  0o_on«-1 

1820  DISP  "Do  you  want  to  (D)raw.  (E)rase,  Plot  with  Pen  (0),  (1),  (2)?" 

1830  LINPUT  Hessagel 

1840  SELECT  Hesaagel 

1850  CASE  "0"."d" 

1860  PLOTTER  IS  3. "INTERNAL” 

1870  GRAPHICS  ON 

1880  PEN  1 

1890  Go  on>1 

1891  CASE  ^0" 

1892  PLOTTER  IS  70S."HPGL” 

1893  GRAPHICS  ON 

1894  PEN  0 

1895  PLOTTER  IS  3. "INTERNAL" 

1896  GRAPHICS  ON 

1898  PEN  1 

1899  Co  on>1 

1900  CASE  ^1" 

1910  PLOTTER  IS  705."HPGL" 

1920  GRAPHICS  ON 

1930  PEN  1 

1940  Go  on>1 

1941  CASE  “2* 

1942  PLOTTER  IS  705. "HPCL" 

1943  GRAPHICS  ON 
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1949  Go  OII>1 

1950  C«SE 

1960  PLOTTER  IS  3.*INTERM«L* 

1970  GRAPHICS  ON 

1980  PEN  -1 

1990  Go  on>1 

2000  CASE  USE 

2010  PRINT  CHR$(7) 

2020  END  SELECT 

2030  UNTIL  Go  on-1 

2040  SUBOID  ~ 

2080  SUB  Got  oil  aoMCElllpao) 

2090  I  ThTa  aubroutlnoa  pronpta  t 

2100  CON  /Elltpsoa/  Natrx(2.2.M0) 


laor  for  tho  Mon  of  an  ollipaa. 
^(20).Xbar(20),Hajor  asia(20), Minor  ax 


la(20).ThoU(20),Anglo(20).Ealn  x(20).EBax  x(20).Eoln_y(20).Ei£x_]r(20) 


2110  Cloaraoroon 

2120  PRINT  TABXX(l,10)tniHAT  IS  X-BAR  for  ELLIPSE  "iElllpao»"7« j 
2130  INPUT  Xbor(Elllpaa) 

2140  PRINT  Xbar(Elllpao) 

2150  PRINT  TABXT(1.11)j*WHAT  IS  I-BAR  for  ELLIPSE  •sElllpaot*7"s 

2160  INPUT  Xbar(Elllpao) 

2170  PRINT  Tbar(Elllpaa) 

2180  Pauaooblt 

2190  SUBEND 

2230  SUB  Got  oooarlanoo(Elllpao) 

2240  I  iftla  routino  gota  tho  oovarianeo  aatrlx  for  an  olllpao 
2250  CON  /Elllpaoa/  Natrx(2,2.40),Xbar(20).Ybar(20),Hajor  axia( 20). Minor  ax 
18(20), Thota(20),Anglo(20).BHln  x(20).teax  x(20).Ealn_]r(20),Eaax_y(20)  ~ 

2260  Cloaraoroon  ~  ” 

2270  FOR  K-1  TO  2 

2280  PRINT  TABXTd ,3*K*6)«*ENTER  ELEMENT  ("!K;","{Ks*)  IN  THE  COVARIANCE 

MATRIX  FOR  ELLIPSE** iElllpao;":  *i 
2290  INPUT  NatrxCK.K.Elllpao) 

2300  PRINT  Natrx(K.K,Clllpao) 

2310  NEXT  K 

2320  PRINT  TABXTd, 3*9) l*ENTER  ELEMENT  (1,2)  IN  THE  COVARIANCE  MATRIX  FOR  E 
LEHENT#"sIj":  "j 

2330  INPUT  Natrxd,2.Elllpao) 

2340  Nntrx(2,1 ,Elllpao)-Natrx(1 .2,Elllpao) 

2350  PRINT  Hatrx(1.2.Elllpao) 

2360  Pauaoablt 

2370  SUBEND 
2410  SUB  Pauaoablt 

2420  I  Pauao  and  wait  for  a  oarriago  rotum 

2430  DISP  "Typo  ENTER  to  eontlnuo...”; 

2440  INPUT  QarbagoS 

2450  SUBEND 
2460  SUB  Cloaraoroon 
2470  I  Cloar  tho  aoroon 

2480  PRINT  CHR$(12) 

2490  SUKND 

2530  SUB  Got  axoa(Elllp8o) 

2540  I  ThTa  aubroutlno  gota  an  oliipao  in  toroa  of  tho  axoa  and  tho 
2550  I  anglo  of  orlontatlon.  Thoao  aro  convortod  Into  a  oovarianeo 

2560  1  aatrlx. 

2570  COM  /Elllpaoa/  Natrx(2.2, 40). Xbar(20),Tbar(20), Major  axla(20).Mlnor_ax 

la(20),Theta(20},Anglo(20),Cain  x(20),Eaax  x(20),Eain  ]r(20),EBiax  y(20)  ~ 

2580  Probjeonatant>2*L0G( )  ~ 

2590  Clovaeroon 

2600  PRINT  "WHAT  IS  THE  RADIUS  OF  THE  MAJOR  AXIS  OF  ELLIPSE  "sElllpaot*  :  " 
t 

2610  INPUT  Major  axla(Elllpao> 

2620  PRINT  Hajor2axla(Elllpao) 

2630  PRINT  *NHAT'~ tS  THE  RADIUS  OF  THE  MINOR  AXIS  OF  ELLIPSE  "lElllpaot"  :  * 


2640  INPUT  Minor juisCElllpM) 

2650  PRINT  Minorjuil3(Ellip8e) 

2660  PRINT  "WHAT~ IS  THE  ANCU  (IN  DECREES)  OF  THE  MAJOR  AXIS  OF  ELLIPSE  "lE 
lllpM;*  :  *; 

2670  INPUT  Antl«(ElllpM) 

2680  PRINT  Ansl«(ElllpM> 

2690  i  Convort  to  Covnrlane*  Matrix 
2700  Theta(Elllpaa)-PI*Ansle(EIXlpaa)/l80 

2710  Matrxd .1  ,EUlpaa)-(Major  axla(Elllpae)*C0S(ThBta(Elllpaa)))‘2*(Mlnor 

axl3(EHlpaa)»SIM(Theta(EUlpaa)))*I  ~ 

2720  Matrx(2,2.EUlp33)-(Major  axl3(Elllpse)*SlN(Th3ta(Elllp33)))*2»(Mlnor 

axls(ElXlp3e)«C0StTh3ta(Elllp3«)))*2 

2730  Matrx(2.1iEXXlpsa)-((Major  axls(EXXlp33))*2-(Mlnor  axl3(EXXlp3e))‘2)*C 

0S(Thota(EXXlp3e))*SIN(Th3U<EXXlp33T) 

2740  Matrxd  ,2.EUlp30)-Matrx(2,1  .EXXlpsa) 

2750  FOR  J-1  TO  2 

2760  FOR  K-1  TO  2 

2770  Matrx( J,K.EXXlpa3)>Matrx(J.K,ElXlp33)/(-Prob  constant) 

2780  NEXT  K  “ 

2790  NEXT  J 

2800  PRINT 

2810  Pauaaablt 

2820  SUBEND 

2860  SUB  DlapXayjeov(EXXlpae) 

2870  I  DlapXaJa  the  covariance  aatrlx  for  an  ellipse 

2880  COM  /Ellipses/  Matrx(2,2. 40), Xbar(20).Ybar(20).MaJor_axls(20), Minor  ax 

ls(20),Thata(20),Angle(20).Enln  x(20).E^x  x(20).Ealn  y(20),Eauc_y(20)  ~ 

2890  rOR  1-1  TO  2  “  “ 

2900  FOR  J-1  TO  2 

2910  PRINT  Matrx(I.J.Elllpse);” 

2920  NEXT  J 

2930  PRINT 

2940  NEXT  I 

2950  SUBEMD 

2990  SUB  A0d_povar lance ( First, Second, SuoBwr) 

3000  1  Tlila  subroutine  can  be  used  to  add  Suoaer-Flrst-Second 

3010  COM  /Ellipses/  Matrx(2,2,40),Xbar(20),Tbar(20),MaJor  ax Is (20), Minor  ax 

ls(20),Theta(20),Angle(20),Ealn  x(20),CBax  x(20),Ealny(20),Eaiax  y(20)  ~ 

3020  INTEGER  I,J  ~  ~ 

3030  FOR  1-1  TO  2 

3040  FOR  J-1  TO  2 

3050  Matrxd, J,SuaMr)-Matrx(I,J,Flrst)*Matrx(I,J, Second) 

3060  NEXT  J 

3070  NEXT  I 

3080  SUBBIO 

3120  SUB  Co«blne_elllpae(Flrst,Second,CoBbo) 

3130  I  This  rwtlne  finds  the  "JPL”  coabinatlon  of  the  First  and  Second 

3140  I  ellipses. 

3150  COM  /Ellipses/  Matrx(2,2, 40), Xbar(20),Xbar(20), Major  axls(20), Minor  ax 

ls(20).Theta(20),Angle(20),Ealn  x(20),Emax  x(20),BBlny(20),EMX  y(20)  ~ 

3160  RCALSx1,Sx2  ~  ~ 

3170  I  Calculate  the  new  covariance  aatrix. 

318O  Invert(Flrst,Flrst*20) 

3190  Invert(Second,Seoond*20) 

3200  Addjeovarlanee( Flrat«20 ,Second*20 ,Coabo-20 ) 

3210  lnvert(CoBbo-20,Coabo) 

3220  I  Calculate  the  new  aean  point. 

3230  Sxl-Nstrxd ,1 ,Flrst«20)*Xbar(Flrat)-Matrx(1 ,2,Flrst*20)*Xbar(Flrst) 

3240  Sxl-Sxl*Matrxd  ,1  ,Second*20)*Xbar(Second)«’Matrx(l  ,2,Seeond*20)*Ybar(Se 

cond) 

3250  Sx2-Matrx( 2 , 1 , First «20 )«Xbar( Flrat) -Matrx( 2 .2, Flrst«20 ) •Ybar (First ) 

3260  Sx2-Sx2*Matrx(2,1 ,Second*20)*Xbar(Second)»Matrx(2.2,Second*20)*Xbar(Se 

cond) 

3270  Xbar ( Coobo ) -Ma tr x( 1,1, Coabo ) *Sxl *Matr x (1,2, Coobo )  *Sx2 

3280  Tbar(Coabo)-Matrx(2,1 ,Coabo)*Sx1«Matrx(2,2,Coabo)*Sx2 

3290  SUBEND 


3330  SUB  DlspjtKtr«Ms(ElllpM) 

3310  i  Tht?  routlm  dtsplaya  th«  axtraa*  x  and  y  valuta  for  tha  apaclflcd 
3350  I  alllpaa 

3360  COH  /Ellipaaa/  Hatrx(2.2.«0),Xbar(20).n>ar(20),HaJor_axia(20).Hlnor  ax 

ta(20),Thata(20),Anilo(20),Ealn  x(20).Emax  x(20),Ealn  y(20)fEMx  y(20)  " 

3370  mm  ■£#  •;Elllpaaj"T  "j  “ 

3360  mm  "Xaln  ";EHln  xCElllpaa); 

3390  PRINT  ■:  Xmax  ";eSuc  xIEUlpaa): 

3900  PRINT  •;  tain  "jEaln^CElllpsa); 

3910  PRINT  "I  laax  "iEaax  y(Elllpae) 

3920  SUBEND 

3930  SUB  Oraw_elllpaa(Elllpae, Xaln. ]Qux. Tain. Xaax.ProbjeonaUnt) 

3990  I  ihiT  routino  drawa  on#  (1)  alllpaa  on  tha  currant  plottar  davice. 
3950  COH  /Ellipaaa/  Matrx(2.2.90).Xbar(20).ybar(20),HaJor  axla(20).Hlnor_ax 

la(20),ThataC20).Angla(20),Ealn  x(20).EBax  x(20),Caln  y(20)iEnix_f(20) 

3960  INTEOER  Sifn  “  ” 

3970  REAL  Xpolnt.Upolnt,Tpoint.Taapl.Taap2,Taap3>Taap9,Taap5 

3980  Gat_bounda(Cllipaa) 

3990  Inv7rt(Elllpaa.Elllpaa*20) 

3500  GRAPHICS  ON 

3510  I  Draw  top  half  of  tha  alllpaa 

3520  Slsn-l 

3530  GOSUB  Orawjialf 

3590  I  Draw  botToa  half  of  tha  alllpaa 

3550  Slgn—I 

3560  GOSUB  Drawjialf 

3570  t  Flnlah  u^ 

3580  HOVE  (Xbar(Elllpaa)-Xaln)*100/(XBax-Xaln),(Tbar(Elllpaa)-Taln)»100/(Xa 

ax-Tain) 

3590  DRAW  (Xbar(Elllpa8)>Xaln)*100/(XBax-XBln).(n>ar(Elllpsa)-yaln)*lOO/(ya 

ax-laln) 

3600  !  Okay,  dokay,  ua'ra  dona  hara. 

3610  SUBEXIT 

3620  Orawjialf:  I  Draw  half  of  tha  alllpaa 

3630  '~i  If  Signal,  than  draw  top  half;  If  Sign— 1,  than  bottom  half. 

3690  FOR  XpolntaEainjc(Elllpaa}  TO  EBaxj((Ellipaa>  STEP  .1 

3650  GOSUB  Draw  1  jpolnt  " 

3660  NEXT  Xpolnt  “ 

3670  XpolnfEBax_x(Elllpaa} 

3680  GOSUB  Draw  l^olnt 

3690  RETURN  “ 

3700  I 

3710  I  Computa  aach  point  and  draw  tha  naw  llna.  (It'a  hara  ainca  wa  call 

3720  I  It  twloa 

3730  1 

3790  0raw_1  joint:  I 

3750  "UpolntaXpoint-XbarlElllpaa) 

3760  Tamp1-Matrx(1,2,Elllpaa*20}*Upoint 

3770  Tomp2>Tampl*Tampt'^trx(2,2,Elllpaa*20)*(Hatrx(1 ,1  ,Elllpaa*20)*Upolnt* 

UpolntaProb  oonatant) 

3780  iT  Taap2<10*(-10}  THEN  Tai92-0 

3790  Xpolnt-(-TaHpi*Slgn*SQR(Tamp2))/Hatrx(2,2,Elllpaa*20)»Xbar(Elllpao) 

3800  IF  Xpolnt<-EHln  x(Elllpaa)a.001  THEN 

3610  HOVE  (Xpolnt>Xain)*100/(Xmax-Xmin},(Xpolnt-Xaln)*100/(Xnax-Xnln) 

3820  ELSE 

3830  DRAW  (Xpolnt-Xmin)*i00/(XwXaln),(Xpoint-Xnin)*100/(Xaax-Xnin) 

3690  END  IF 

3650  RETURN 

3660  SUBEND 

3900  SUB  Taat(Flrat,Saoond,Mork,Taat_val} 

3910  i  Thla  routlnaa  oaloulataa  tha  accaptanca  taat  erltarla  for  Flrat  and 

3920  I  Saoond.  Hork  la  uaad  aa  a  "scratchpad". 

3930  COH  /Ellipaaa/  Natrx(2,2,90),Xbar(20),Xbar(20).Hajor_axia(20),Nlnor  as 

ls(20),Thata(20),Angia(20),Enlnjc(20).Enax  x(20).Enin jr(20>,Emaxjr(20)  -  " 

3990  Add jovarianoaC  First,  Saoond, WorkT 

3950  lnv7rt(Mork,WorH«20) 
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